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Measurements  of  a  1230-km  baseline  were  made  during  an  eight-week  period  in  the 
fall  of  1987  using  Trimble  *l(JUOSX  single-frequency,  five  channel  Global  Positioning 
System  (GPS)  receivers.  Twenty-eight  days  of  carrier  phase  data  were  processed  using 
correlated  triple  differences  with  fixed  satellite  orbits,  the  broadcast  ephemerides,  a 
modified  Hoplield  tropospheric  model,  and  without  ionospheric  correction  to  determine 
the  accuracies  and  precisions  of  the  slope  distance  and  baseline  components.  The  data 
were  processed  in  ever  increasing  observing  sessions  to  determine  the  optimized  obser¬ 
vation  periods  required  to  achieve  various  orders  of  geodetic  accuracies. 

The  accuracies  of  the  slope  distances  were  better  than  1.0  ppm  for  any  observing 
period.  The  day-to-day  repeatabilities  of  the  slope  distance  measurements  were  better 
than  1.0  ppm  (2cr>  after  20  minutes  of  observations.  Accuracies  and  repeatabilities  (2<j) 
of  the  baseline  components  were  better  than  10.0  ppm  after  20  minutes  of  observations. 
The  correlated  triple  difference  results  were  on  the  order  of  previous  GPS  surveys  that 
used  higher  resolution  differencing  or  external  timing  aids.  Discussions  include  the  ef¬ 
fects  of  ephemeris,  tropospheric  and  ionospheric  errors,  and  dilution  of  precision. 

Observation  periods  and  mean  slope  distance  errors  were  reduced  when  observations 
started  close  to  and  included  the  infinite  peak  of  the  Position  Dilution  of  Precision 
(PDOP).  The  smallest  variances  were  associated  with  observations  about  the  infinite 
PDOP  peak. 
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I.  INTRODUCTION 


.9  priori  knowledge  of  the  observation  periods  required  to  achieve  specified  orders 
of  geodetic  accuracies  is  important  in  planning  efficient  and  productive  geodetic  surveys. 
While  terrestrial  surveys  have  specified  field  and  processing  procedures  and  standards  to 
categorize  the  geodetic  accuracies  of  surveys  [Federal  Geodetic  Control  Committee. 
19S4),  only  recently  have  standards  been  proposed  for  surveys  conducted  with  the  Global 
Positioning  System  (GPS)  [Federal  Geodetic  Control  Committee.  1986].  Among  the 
proposed  requirements  are  standards  for  the  length  of  observing  periods  and  satellite 
geometry. 

Field  studies  by  Remondi  [1984]  and  numerical  simulations  by  Fell  [19 SO | ,  Langley 
et  al.  [1984],  and  Landau  and  Eissfeller  [1986]  studied  optimized  observation  periods,  but 
for  baselines  less  than  100  km.  Cannon  et  al.  [  1985],  Bock  et  al.  [1984],  Goad  et  al. 
[1985],  Mader  and  Abell  [1985],  and  Bertiger  and  Lichten  [1987]  conducted  long  baseline 
surveys,  but  did  not  study  optimized  observation  periods.  One  of  the  objectives  of  this 
thesis  is  to  fill  the  gap  between  the  above  studies,  /.<?.,  examine  the  optimized  observatior 
periods  for  a  long  baseline. 

The  optimized  times  will  be  examined  using  the  correlated  triple  difference  carrier 
beat  phase  observable  because  of  its  insensitivity  to  integer  ambiguities  and  loss  of  lock 
of  the  GPS  carrier  by  the  receiver.  Another  objective  of  this  thesis  is  to  add  to  the  body 
of  triple  difference  accuracy  testing  following  a  recommendation  by  Remondi  [1984.  p. 
259J:  "More  testing  is  required  to  establish  the  full  accuracy  potential  of  the  triple  dif¬ 
ference  method. 

GPS  carrier  phase  and  pseudorange  measurements  were  made  during  an  eight-week 
period  in  the  fall  of  1987  with  Trimble  4000SX  single-frequency,  5  channel  receivers. 
The  long  baseline  is  approximately  1230  km  in  length  between  the  National  Oceanic  and 
Atmospheric  Administration's  (NOAA)  Western  Regional  Center  located  at  Sand  Point 
in  Seattle,  Washington,  and  the  Naval  Postgraduate  School  (NPS).  Monterey, 
California.  The  baseline  was  determined  by  locating  the  positions  of  its  ends  by  con¬ 
necting  them  by  independent  short  baseline  surveys  from  nearby  Very  Long  Baseline 
Interferometry  c  V L B I )  horizontal  control  points.  The  results  of  those  surveys  form  the 
reference  to  which  accuracy  will  be  determined. 


Additionally,  studies  for  repeatability  were  conducted  following  another  recommen¬ 
dation  by  Remondi  (1984,  p.  263]  to  enhance  the  capabilities  of  GPS  measurements. 
The  recommendation  was  to  perform  extensive  repeatability  studies  on  non-vanme 
baselines  for  verifying  and  improving  the  GPS  modelling. 


II.  BASELINE  DETERMINATION  USING  GPS 


A.  INTRODUCTION 

The  Department  of  Defense  s  <DOD)  Global  Positioning  Svstem  is  intended  to 
provide  accurate  positioning  and  precise  tuning  for  navigation  purposes  by  broadcasting 
codes  superimposed  on  two  radio  carrier  frequencies  from  satellites.  The  satellites  are 
placed  in  a  constellation  so  that  at  least  four  satellites  arc  visible  globally.  The  Precise 
Code  (P  code)  will  be  limited  to  authorized  DOD  users.  The  Coarse  Acquisition  (C  A) 
code  provides  real-time  accuracies  to  about  100  m  [Baker,  19S6]  and  is  available  to 
anyone. 

The  codes  provide  their  transmit  times,  satellite  orbit  and  clock  information,  and 
information  to  enable  any  receiver  to  lock  onto  other  GPS  satellites.  The  orbital  infor¬ 
mation  (ephemeris)  provides  the  position  of  the  satellite.  The  receiver  measures  the  time 
delay  between  the  receipt  of  the  C  A  code  and  its  transmission  time.  The  time  delay  can 
be  transformed  into  an  apparent  slant  range  from  the  satellite  s  known  position  to  de¬ 
termine  the  location  of  the  receiver.  Since  it  includes  delays  due  to  receiver  clock  errors 
and  the  effects  of  atmospheric  refraction,  the  apparent  slant  range  is  referred  to  as  the 
pu'uJorange.  A  minimum  of  four  satellites  are  required  to  solve  the  system  of  range 
equations  for  the  receiver's  coordinates  and  clock  errors. 

While  the  C  A  code  provides  real-time  location,  it  does  not  meet  the  accuracy  re¬ 
quired  for  precise  geodetic  work.  Nevertheless,  GPS  makes  possible  a  higher  resolution 
via  the  carrier  signal.  Though  the  carrier  itself  does  not  contain  the  orbital  and  timing 
information,  which  would  have  to  be  supplied  by  some  other  means,  it  does  offer  a 
higher  resolution  because  of  its  19-cm  wavelength 

B.  THE  MONTEREY-SAND  POINT  BASELINE 
1.  General 

The  length  of  the  Monterey-Sand  Point  baseline  (Figure  1)  was  computed  by 
ddferencing  the  World  Geodetic  System  1984  (WGS84)  [Defense  Mapping  Agency, 
1987]  Cartesian  coordinates  determined  for  Monterey  and  Sand  Point  by  short  baseline 
GPS  surveys  from  known  horizontal  control  points.  The  precision  and  agreement  with 
terrestrial  survey  results  of  short  baseline  GPS  measurements  using  single  frequency, 
double  difference  solutions  are  well  documented  (eg.,  [Remondi,  1984],  [Goad  and 
Rcmondi,  1984|,  [Bock  et  al.,  1984]). 
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Figure  1.  Monterey-Sand  Point  Baseline  and  environs:  Insets  not  drawn  to  scale. 

At  Sand  Point,  on-site  meteorological  measurements  were  made  near  the  middle 
of  the  observing  session.  For  the  Monterey  antenna  determination  meteorological 
measurements  were  made  every  half  hour  and  the  mean  of  all  the  measurements  was 
used  in  the  processing.  The  Trim640  solutions  were  obtained  using  uncorrclated  double 
differences  and  estimating  initial  integer  ambiguities.  The  ambiguities  were  fixed  to  the 
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integer  values  that  produced  the  smallest  residuals.  A  tropospheric  factor  was  estimated 
along  with  the  integer  ambiguities  and  the  baseline  components  in  the  least-squares 
processing. 

The  horizontal  control  points  used  for  the  reference  stations  in  determining  the 
coordinates  of  the  antennas  were  mobile  Very  Long  Baseline  Interferometry  (VLBI) 
sites.  The  NADS3  Cartesian  coordinates  for  the  YLB1  sites  were  provided  by  the  the 
Gravity,  Astronomy  and  Space  Geodesy  Branch  of  the  National  Geodetic  Survey  tNGSi 
[Abell.  19S7],  The  NADS3  coordinates  were  determined  in  August  19S7  from  a  global 
adjustment  of  VLBI  surveys.  Carter  et  al.  [1985)  described  the  NGS  VLB1  program. 
The  Defense  Mapping  Agency  Hydrographic  Topographic  Ollice  (DMAIITCi  validated 
the  direct  transformation  of  the  VLBI  Cartesian  coordinates  to  WGSS4  Cartesian  coor¬ 
dinates  [Kumar,  19S8|. 

2.  Monterey  Coordinates 

The  coordinates  for  the  Monterey  antenna  location  were  determined  by  aver¬ 
aging  two  surveys  conducted  on  separate  days  from  the  VLBI  site  FT  ORD  NC.MN 
19S1.  Table  1  lists  the  observing  sessions  used  to  determine  the  WGSS4  coordinates  of 
the  Monterev  antenna.  Table  2  lists  the  results  of  the  GPS  survevs. 


Table  I.  MONTEREY  ANTENNA  LOCATION  SURVEYS:  From  FT  ORD 
NCMX  1981.  PST  =  Pacific  Standard  Time. 


Date 

Start 

lime 

(PST) 

End 

Time 

(PST) 

Number  of 
Double 
Dilference 
Observations 

Percent 

Rejected 

Slope 

Distance  (m) 

°b 

09  16  S7 

10  IS 

1210 

996 

7 

12139.8125 

0.025 

09  18  87 

0914 

1120 

1136 

4 

12139.8119 

0.025 

The  Monterey  coordinates  were  computed  from: 
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Monterey  coordinate 
Ft  Ord  coordinate 

Variance-weighted  mean  baseline  component 
The  uncertainties  of  the  Monterey  coordinates,  <tc  ,  were  computed  from: 
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Table  2. 


DATE 


ID  16  S'1 
(D  IS  ST 


RESULTS  OF  MONTEREY  ANTENNA 
SURVEYS:  WGS84  Cartesian  Coordinates  (meters). 


LOCATION 


FT  ORD  X 


li 


697026.493 
09  '026.493 


aORD  x 


0.007 
0.00  7 


MRY  X 


-2707340.O93 

-2707340.095 


DATE 

FT  ORD  Y 

aORD  Y 

AY 

09  16  87 
09  18  87 

-4354393.309 

-4354393.309 

0.010 

0.010 

917.693 

917.689 

DATE 

FT  ORD  Z 

a  ORD  Z 

AZ 

09  16  87 
09  IS  S7 

37S8077.77S 

3~SS077.778 

0.009 

0.009 

-6337.391 

-6337.389 

Mean:  -27O~340.o94 


MRY  Y 


-4353475.617 

-4353475.620 


Mean:  -4353475.61S 


MRYZ 


0.041  3781740. 3S7 

0.044  37S  1740.389 


Mean:  3781740.388 


-> 
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where  <Jr0ftD  is  the  uncertainty  of  the  Ft  Ord  coordinate. 

One  month  prior  to  the  surveys  originating  from  Ft  Ord,  two  other  surveys  were 
conducted  from  the  satellite  Doppler  horizontal  control  point  NAVAL  POST  GRAD 
31965  DOPPLER.  The  Doppler  station  is  approximately  300  m  north  of  the  Monterey 
antenna  location  (Figure  1).  The  double  differenced  GPS  carrier  phase  solutions  of  the 
two  Doppler-originating  surveys  yielded  the  following  Monterey  coordinates: 

X  -2707339.725  m 

Y  -4353475.654  m 

Z  3781740.264  m 

The  baseline  components  had  an  uncertainty  of  ±0.002  m.  but  the  Doppler  station  has 
an  uncertainty  of  ±2  m  in  each  coordinate  before  transformation  to  WGS84.  The  results 
of  the  Doppler  surveys  were  not  used  in  determining  the  Monterey  coordinates  because 
of  the  large  uncertainty  in  the  Doppler  station  location.  The  three-dimensional  positions 


of  the  Monterey  antenna  from  the  Doppler  and  the  VLBI  originating  survey  agree  to 
better  than  0.5  m. 


The  two  days  of  pseudorange  observations  at  the  Montere}  receiver  were  each 
subjected  to  a  least-squares  estimation  and  then  averaged  together  to  yield  the  Monterey 
coordinates: 

X  -2707334.248  ±1.3  m 

Y  -4353466. Su7  ±0.6  m 

Z  3781741.330  ±2.2  m 

where  the  two-day  standard  deviations  are  given.  The  deviation  of  the  pseudorange 
from  the  differenced  carrier  phase  derived  Monterey  coordinates  is  expected  because  of 
the  coarser  resolution  of  the  C  A  code  and  because  the  pseudoranges  are  corrected  for 
neither  tropospheric  nor  ionospheric  delays. 

3.  Sand  Point  Coordinates 

The  Sand  Point  antenna  coordinates  were  determined  by  one  90-minute  GPS 
survey  (Table  3)  from  the  mobile  VLBI  site  Aviation  2  which  is  530  m  distant 
(Figure  1). 


Table  3.  RESULTS  OF  SAND  POINT  ANTENNA  LOCATION 
SURVEY:  From  double  difference  carrier  phase  solutions.  Distances 
and  WGS84  Cartesian  Coordinates  are  in  meters. 


Aviation  2  Coordinate 

&  in 

A 

U 

Sand  Point  Coordinate 

X  =  -2295347.760 

-408.361 

Bj|S0 

-2295756.121 

0.017 

Y=  -3638029.429 

MtHJMP; 

ragfifB 

-3637699.228 

0.028 

Z=  4693408.964 

ihjoeem 

73.813 

0.003 

4693482.777 

0.032 

The  Sand  Point  antenna  coordinates  computed  from  the  pseudorange  data 

were: 

X  -2295749.623  m 

Y  -3637694.279  m 

Z  4693488.472  m 

and  standard  deviations  are  not  given  because  this  is  a  single  observing  session  and  the 
program  does  not  provide  solution  standard  deviations. 


4.  Monterey-Sand  Point  Baseline  Components 

First  the  Monterey-Sand  Point  baseline  components  were  computed  by  sub¬ 
tracting  the  Sand  Point  Cartesian  coordinate  from  the  respective  Monterey  Cartesian 
coordinate.  The  uncertainties  in  the  baseline  components  were  computed  as  the  square 
root  of  the  sum  of  the  squares  of  the  uncertainty  of  the  Monterey  and  Sand  Point  co¬ 
ordinates.  The  slope  distance  was  computed  as  the  square  root  of  the  sum  of  the  squares 
of  the  baseline  components.  The  slope  distance  uncertainty  was  computed  as: 

A  A’  2  ,  ,  A  Y  ,2  ,  ,  AZ  J 

°b  =  (  —  Vsx)  +(  —  >  +(  —  ej±z) 


•  2  ,  ,  AT  .2  .  AZ  2 

"ax)  +(  —  +1  —  °AZ> 


where  b  is  the  slope  distance.  Using  the  information  from  Tables  2  and  3  gives 
Table  4,  which  is  used  as  a  standard  to  estimate  accuracy. 

Table  4.  STANDARD  BASELINE  DISTANCES  AND  2-SIGMA  VALUES. 


AX 

-411583.973 

+  0.080  m 

(0.19  ppm) 

AY 

-715776.390 

±  0.116  m 

(0. 16  ppm) 

AZ 

-911742.388 

±  0. 109  m 

(0.12  ppm) 

Slope  Distance 

1230045.280 

±  0.109  m 

(0.09  ppm) 

C.  THE  ONE  WAY  CARRIER  BEAT  PHASE 

The  development  of  the  carrier  beat  phase  technique  and  a  model  for  its  application 
are  given  by  Remondi  (1984|.  The  phase  measurement  is  done  by  beating  the  received 
carrier  with  a  local  oscillator  internal  to  the  GPS  receiver.  The  slant  range  from  a  GPS 
receiver  to  a  GPS  satellite  can  be  modelled  in  terms  of  the  time  it  takes  the  signal  to 
travel  or  the  number  of  cycles  that  occur  between  the  satellite  and  the  receiver.  The 
range  in  cycles  will  consist  of  an  integer  and  fractional  number  of  cycles.  When  a  GPS 
receiver  locks  onto  the  carrier  signal,  it  can  immediately  measure  the  fractional  part  and 
begin  counting  subsequent  integer  cycles,  but  it  cannot  measure  or  account  for  the  initial 
integer  number  of  cycles  that  preceded  the  initial  fractional  part.  These  missing  cycles 
bias  subsequent  measurements  and  are  called  the  initial  integer  ambiguity  biases. 

The  signal  does  not  take  a  direct  path  to  the  receiver  as  it  is  refracted  by  the 
ionosphere  and  troposphere.  Additional  errors  are  caused  by  the  satellite  deviating  from 
its  predicted  orbit,  errors  in  the  satellite  clock,  and  errors  in  the  receiver  clock.  Anti¬ 
cipating  that  the  observables  will  be  used  in  relative  positioning,  that  a  single  frequency 
receiver  will  be  used  and  ignoring  other  error  sources,  such  as  multipath,  the  one-way 
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carrier  beat  phase,  d>b(r.s.i),  observed  at  receiver,  r,  from  satellite.  5,  at  observation 
epoch,  t,  can  be  modelled  to  first-order  as: 

<j>b(r,s,i)  =  (ps(i)  -  4>r(r,i)  -  ~  p(r,s,i ) 

+  CMC*  +fr  -  T  P{r'5'‘)l  ' "  ! ' 

-  r(r./')  +  .-ffr.s.l) 

where: 

i  Epoch  identifier 

r  Receiver  identifier 

s  Satellite  identifier 

Phase  of  received  carrier  signal 
<t>\r,i)  Phase  of  receiver  generated  carrier  signal 
f  Transmitted  frequency  of  carrier  signal 

p(r.s,i)  Satellite-receiver  slant  range 

p(r,s.i)  Time  rate  of  change  of  slant  range 

c  Speed  of  light 

f  Receiver  generated  carrier  frequency 

r(r,i)  Tropospheric  delay 

.•f(r,s,l)  Initial  integer  ambiguity 

...  ...  j(t) 

Sl  EO  <5  (0  T 

l. 

C(2  .o-M  +  ^Sr- 

i. 

c(i)  Mean  clock  offset  for  both  receivers 

*5(/)  Clock  drift  between  both  receivers 

and  parentheses  do  not  indicate  factors  or  functions,  but  simply  enclose  identifiers. 
Brackets  indicate  factors. 

D.  DIFFERENCING  THE  ONE-WAY  CARRIER  BEAT  PHASE 
1.  Single  Difference 

The  single  difference  ( SD )  is  formed  by  differencing  carrier  beat  phase  observa¬ 
bles  from  two  receivers  at  the  same  observation  epochs.  Following  Remondi  [  1984], 
taking  the  difference,  expanding  the  C(r.i')  terms  of  Equation  (2.1),  and  expressing  dif¬ 
ferences  in  the  between-satellite  and  between-receiver  phases  as  fS(i)  gives: 
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SD(s,i )  =  <J>t{2.s.i)  -  (pb{\,s.i) 

=  —■  C/H2.5.I)  -  p(l.s./)] 

+  7-  li>[2.s.i)  -  p([.5./)]c(/) 

( 2.2) 

-  C/i  +/A‘) 

-I - 7"  Cp(2,5,/)  +  1 .5, /)]<)(/) 

+y;[r(2,/)  -  T(l./)]  +  ,f(2.s,l)  -  ,f(l.s.l) 

where  the  terms  are  the  same  as  Equation  (2.1).  Single  differencing  reduces  or  eliminates 
satellite  orbital  and  clock  errors  because  the}’  are  common  to  both  receivers. 

2.  Double  Difference 

A  double  difference  (DD)  is  formed  by  differencing  single  differences  between  a 
reference  satellite,  h,  and  another  satellite  at  the  same  epoch: 

DD(h,s,i)  =  SD{h,i)  -  SD(s,i) 

Because  the  differences  at  the  same  epoch  are  taken  with  the  same  reference 
satellite,  the  double  differences  for  each  epoch  are  correlated.  The  advantage  of  the 
double  differences  is  that  the  clock  dependent  terms  -  ft6(i)  and  [f2  -/]£(/)  *  are  elim¬ 
inated.  The  significance  of  the  removal  of  those  terms  is  to  reduce  from  nanoseconds 
to  microseconds  the  timing  accuracy  required  to  achieve  one  cycle  accuracy.  The 
Trimble  4000SX  achieves  sub-nucrosecond  accuracy  by  using  the  C  A  code  timing  in¬ 
formation  [Ashjaee,  1985|. 

3,  Triple  Difference 

A  triple  difference  (TD)  is  formed  by  differencing  the  double  differences  for  the 
same  satellite  pair  at  some  integer  number  of  succeeding  epochs,  l : 

TD{{)  =  DD(h,s,i  +  /)  -  DD{h,  s,i) 

The  advantage  of  the  triple  difference  is  that  it  eliminates  all  the  time  inde¬ 
pendent  terms,  namely  the  initial  integer  ambiguities,  /f(r,s,I),  and  becomes  insensitive 
to  the  initial  ambiguities  and  any  cycle  slips  when  the  receiver  loses  lock. 

The  disadvantages  of  the  triple  difference  are:  another  level  of  correlation,  loss 
of  resolution,  and  reduced  number  of  observations.  The  triple  differences  are  already 
correlated  with  respect  to  satellite  because  of  the  underlying  double  differences,  and  are 
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further  correlated  with  respect  to  time  because  consecutive  triple  dilTerenced  observa¬ 
tions  will  have  the  DD(h,s,i  +  f)  term  in  common. 

For  short  baselines,  integer  ambiguities  can  easily  be  resohed  because  unmod- 
elied  errors  are  highly  correlated  between  the  two  antenna  sites  and  are  mostly  elimi¬ 
nated  by  the  differencing.  Algorithms  can  take  advantage  of  the  integer  nature  of  the 
initial  ambiguities  and  solve  for  them.  At  longer  baselines  the  unmodelled  errors  are  not 
as  highly  correlated  and  not  eliminated  by  the  dill'erencing.  lhese  errors  fold  into  the 
initial  ambiguities,  so  that  the  ambiguities  are  no  longer  integers.  In  some  cases  the 
ambiguities  cannot  be  resolved  ( e.g .,  Henson  and  Collier,  1986,  Tables  1  and  2). 

Because  of  the  advantages  the  triple  difference  offers  for  long  baselines,  I  use  the 
triple  difference  scheme.  While  the  triple  difference  can  be  decorrelated  by  forming  a 
weight  matrix  [Remondi,  1984],  only  the  correlated  triple  difference  software  was  avail¬ 
able  to  me. 

E.  ERROR  EFFECTS 

For  long  baseline  GPS  surveys,  the  primary  errors  are  satellite  orbit  errors, 
ionospheric  and  tropospheric  delays  [Remondi,  1984,  and  Beutler  et  al.,  1986].  Orbit 
(ephemeris)  errors  are  the  result  of  the  departure  of  the  satellite  from  the  broadcast 
ephemeris  orbit.  The  ephemeris  is  a  predicted  orbit  for  the  satellite.  Orbit  errors  prop¬ 
agate  directly  into  the  baseline  measurements  when  the  GPS  orbit  coordinates  are  fixed 
in  the  processing  [Hothem  and  Williams,  19S5].  Orbit  errors  can  be  the  dominant  error 
source  affecting  the  repeatability  of  long  baseline  measurements  [Lichten  and  Border, 
1987], 

The  magnitudes  of  baseline  errors  increase  with  increasing  baseline  length  because 
of  the  increasing  projection  of  the  ephemeris  error  onto  the  baseline  component  [Fell, 
19S0].  Estimates  of  the  broadcast  ephemeris  error  range  from  25  m  [Beutler  et  al.,  19S6] 
to  100  nr  [Wells  et  al.,  198o].  The  magnitude  of  the  effect  of  the  ephemeris  error  on 
baseline  accuracy  has  been  traditionally  approximated  [Beutler  et  al.,  1984,  Equation 
(2.1)]: 


where: 

b  baseline  length 

p  slant  range  (receiver  to  satellite) 


Eb  error  in  baseline  length 

Ee  error  in  slant  range 

The  slant  range  to  a  GPS  satellite  is  about  20,000  km.  which  translates  to  a  baseline  er¬ 
ror  ranging  from  1  ppm  to  4  ppm  using  Equation  (2.3). 

It  is  expected  that  the  ephemeris  error  for  the  Monterev-Sand  Point  baseline  will  be 
towards  the  lower  end  of  the  range  because  the  ephemendes  are  uploaded  prior  to  the 
satellites  entering  the  Yuma  Proving  Ground  [Russell  and  Schaibiy,  1980J  near  the 
California  border.  The  ephemeris  linearization  error  specification  is  to  1  m  per  day 
[Wells  et  at .,  19$6|. 

Hie  ionosphere  disperses  the  code  (the  group  velocity)  from  the  carrier  phase  (phase 
velocity)  because  the  C  A  code  has  a  frequency  of  1  MHz  while  the  carrier  signal  upon 
which  it  is  superimposed  has  a  frequency  of  1575  MHz.  The  effect  is  to  increase  the 
pseudoranges,  but  decrease  the  carrier  phase  derived  ranges  [Smith,  1987],  Field  exper¬ 
iments  [Beutler  et  at ..  1986|  showed  that  ionospheric  dispersion  shortens  baselines  on  the 
order  of  a  few  tenths  to  perhaps  2  ppm. 

Ionospheric  error  is  proportional  to  the  Total  Electron  Content  (TEC)  along  the 
signal  path  and  the  cosecant  of  the  elevation  angle  of  the  satellite  [Snuth,  1987],  Thus 
the  error  is  greatest  for  low  elevation  angles  and  least  at  the  zenith.  Wells  et  at.  [1986] 
estimated  that  the  range  errors  due  to  the  ionosphere  are  from  150  m  at  the  horizon  to 
5».)  ni  at  the  zenith. 

Ionospheric  activity  is  a  function  of  latitude,  longitude,  time  of  day,  season,  and 
sunspot  activity.  Ionospheric  activity  increases  towards  the  equator  and  towards  the 
sunlit  portions  of  the  earth.  Diurnally,  it  has  a  minimum  near  0600  local  time  with  a 
maximum  around  1600  local.  Ionospheric  activity  increases  with  the  peak  in  the  sunspot 
cycle.  The  minimum  in  the  current  1 1-year  sunspot  cycle  occurred  in  19S6.  Upon  these 
systematic  characteristics  sporadic  ionospheric  disturbances  are  superimposed. [Henson 
and  Collier,  1986] 

The  tropospheric  error  is  proportional  to  the  refractivity  along  the  satellite-receiver 
path  and  proportional  to  the  cosecant  of  the  elevation  angle  [Martin,  19S0[.  The  index 
of  refraction,  n,  is  the  ratio  of  the  speed  of  light  in  a  vacuum  to  the  speed  of  light  in  a 
particular  medium,  in  this  case  the  troposphere.  Because  the  refractive  index  is  a  small 
fraction  greater  than  1.0,  a  more  convenient  unit  to  work  with  is  refractivity,  .V,  where 
.V  =  (/i  -  1)  x  106  .  The  magnitude  of  the  tropospheric  biases  range  from  20  m  for  10° 
elevation  angles  to  2.3  m  at  the  zenith  [Wells  et  at.,  1986). 
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Henson  and  Collier  (1986]  have  shown  that  triple  difference  measurements  are  un¬ 
affected  by  path-dependent  ionospheric  bias  errors,  but  by  path-dependent  Gradient 
ionospheric  errors,  and  \Iartin  (19S0|  has  estimated  that  the  combined  ionospheric  and 
tropospheric  gradient  errors  are  on  the  order  of  meters  per  hour  and  are  proportional 
to  the  cosecant  of  the  elevation  ancle. 
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hi.  data  collection,  processing,  and  analysis 


a.  EQUIPMENT  CONFIGURATION 

1.  Hardware 

A  complete  description  of  the  Trimble  -iuooSX  receiver  is  given  by  Trimble 
Navigation  (1987a).  NPS  operates  three  Trimble  4uuoSX  GPS  Survey  or  receivers,  of 
which  two  were  used  in  this  study.  The  40noSX  is  capable  of  observing  the  C  A  code, 
integrated  Doppler  and  carrier  beat  phases  of  up  to  five  satellites  simuitaneousiv.  Its 
ability  to  use  the  C  A  code  allows  the  receiver  to  be  used  as  a  stand  alone  navigation 
system  which  determines  position  using  Doppler-smoothed  pseudoranges  and  velocities 
[Ashjaee,  1985). 

For  precise  relative  positioning  the  400USX  can  transmit  its  data  through  an 
RS232  port  to  a  microcomputer  for  storage  on  Poppy  disk,  for  post-processing.  The 
dwOSX's  ability  to  use  the  C  A  code  allows  it  to  decode  the  GPS  navigation  messages 
so  that  it  can  track  satellites  automatically  and  determine.  Most  importantly,  it  uses  the 
C  A  code  in  a  time  transfer  mode  to  determine  any  olTset  and  drift  of  its  own  clock  and 
thus  provide  accurate  time  tags  for  the  observations  without  an  external  atomic  clock 
or  synchronization  with  the  receiver  at  the  other  end  of  the  baseline. 

The  receivers  were  left  on  continuously  to  allow  unattended  data  collection. 

Multipath-resistant  Trimble  microstrip  antennas  were  installed  at  both  the 
Monterey  and  Sand  Point  locations. 

2.  Software 

For  relative  positioning,  the  receiver  is  controlled  from  the  microcomputer  by 
Version  D  of  Trimble  s  Datalogger  program.  The  reference  position  (the  gcoderic  coor¬ 
dinates  of  the  antenna)  and  the  particular  options  chosen  must  be  entered  into  the  re¬ 
ceiver  via  the  receiver  keypad.  The  receivers  were  set  to  use  the  reference  position  height 
for  point  positioning  when  less  than  four  satellites  were  available. 

Each  observation  session  was  initialized  to  log  data  when  a  rrummutn  of  four 
satellites  were  15°  above  the  antenna's  horizon.  Five  satellites  were  designated  for  each 
observing  session.  The  software  logs  the  observables  and  receiver  clock  parameters  to 
a  floppy  disk  every  15  seconds  and  the  C  A  code-determined  antenna  position  and  Po¬ 
sition  Dilution  of  Precision  (PDOP)  every  five  minutes.  The  GPS  navigation  message 
is  logged  to  a  separate  file  at  the  beginning  of  the  session. 


B.  SATELLITE  OBSERVATION  PLAN 

1.  Satellite  Selection 

The  same  live  satellite  (SVl  <(■>.  9.  11,  12.  and  13)  were  used  for  the  entire  cght- 
week  observation  period.  These  live  satellites  were  visible  at  both  stations  l'or  over  l» »U 
nunutes,  and  of  these  live  satellites  four  were  visible  for  three  hours  (Figure  2). 

START  DATE/TIME:  1987/10/17  15:21:15.  DAY  OF  YEAR  290  TOW  573675. 

STOP  DATE/TIME:  1987/10/17  18:16:  1.  DAY  OF  YEAR  290  TOH  584160. 

DATA  AVAILABLE: 

STATION:  SAN0  POINT 

SV  6| .  | 

SV  9  | . | 

SV  11  | . | 

SV  12  | . | 

SV  13  I . | 

STATION:  MONTEREY 

SV  6  I .  | 

SV  9  | . | 

SV  11  | . | 

SV  12  |  1 

SV  13  | . | 

Figure  2.  Satellite  availability:  Each  dot  represents  10  observations;  each  col¬ 
umn,  10  epochs.  From  Trim640  output. 

2.  Position  Dilution  of  Precision  (PDOP) 

Trimble  Navigation  [1987b]  recommends  that  observations  include  the  time  that 
the  PDOP  goes  to  infinity.  The  Federal  Geodetic  Control  Committee  [1986]  notes  that 
initial  results  from  investigations  indicate  that  best  results  may  be  achieved  when  the 
Geometric  Dilution  Of  Precision  (GDOP)  is  changing  value  during  the  observing  ses¬ 
sion.  and  proposes  that  observing  sessions  start  with  a  high  GDOP  and  stop  with  a  low 
GDOP. 

Position  Dilution  of  Precision  is  a  component  of  the  GDOP.  GDOP  is  a 
measure  of  how  satellite  geometry  degrades  point  position  accuracy  [Jorgensen,  1984]. 
For  computational  ease  in  the  navigation  solution  the  GDOP  is  defined  as  the  square 
root  of  the  trace  of  the  covariance  matrix  of  the  errors  in  position  and  time  with  the  range 
errors  set  to  one  [Milliken  and  Zoller,  1980].  The  role  and  definition  of  GDOP  and 
PDOP  in  GPS  point  positioning  were  applied  to  GPS  relative  positioning,  i.e.,  good 


PDOP  would  provide  better  accuracy  than  poor  PDOP  [King  ei  al.  1985).  Good  and 
poor  PDOP  are  shown  in  Figure  3. 

Landau  and  Eissfeller  [1986|.  using  numerical  simulations  in  which  they  assumed 
a  full  IS  satellite  constellation  and  a  receiver  that  could  track  those  satellites  that  mini¬ 
mized  GDOP,  found  that  better  accuracy  lor  triple  difference  solution  lor  a  68  km 
baseline  generally  corresponded  to  high  GDOP.  They  used  a  more  complete  GDOP  that 
included  consideration  of  ionospheric,  tropospheric  and  satellite  position  errors  which 
are  neglected  in  the  conventional  GDOP. 

The  40UUSX  does  not  record  GDOP,  but  it  does  record  PDOP  every  five  min¬ 
utes.  PDOP  relates  to  GDOP  as:  GDOP2  =  PDOP2  +  TDOP2,  where  TDOP  is  the  Time 
Dilution  Of  Precision,  the  error  in  the  user  clock  bias  multiplied  by  the  velocity  of  light. 
The  expected  uncertainty  in  a  GPS  point  positioning  solution  is  a  product  of  the  PDOP 
and  the  expected  slant  range  error.  The  difference  between  PDOP's  at  Monterey  and 
Sand  Point  remained  less  than  1.0  for  the  entire  eight  week  observing  period.  The  PDOP 
peaks  at  two  times  in  an  observing  session  (Figure  4)  -  60  minutes  and  150  minutes. 
The  PDOP  peaks  occur  near  when  the  satellites  lie  in  a  common  plane  causing  the  sol¬ 
ution  of  the  linearized  range  equations  to  diverge  and  the  PDOP  becomes  infinite 
[Jorgensen,  1984). 

C.  METEOROLOGICAL  PARAMETERS 

Meteorological  parameters  are  needed  for  the  tropospheric  correction  model  used  in 
the  processing  software.  On-site  meteorological  observations  were  not  available  for 
Sand  Point.  Instead  observations  from  the  Weather  Service  Office  (WSO)  of  the 
Seattle-Tacoma  International  Airport  (SEA-TAC)  were  used  for  Sand  Point.  SEA-TAC 
is  approximately  27  km  from  Sand  Point  (Figure  1).  The  NPS  Department  of  Meteor¬ 
ology  routinely  collects  real-time  hourly  observations  of  sea-level  pressure,  temperature 
and  relative  humidity  from  the  National  Weather  Service's  data  network.  Observations 
that  pertained  to  GPS  observing  sessions  were  entered  into  a  file  that  is  accessed  by  the 
batch  file  building  program. 

While  the  Federal  Geodetic  Control  Committee  (1986)  has  proposed  the  use  of  on¬ 
site  meteorological  parameters,  researchers  have  had  success  using  standard  atmosphere 
parameters  for  satellite  geodesy  [Fell,  1976]  or  extrapolated  meteorological  data 
[Rathacher  et  al.,  1986], 
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Figure  3.  Poor  and  good  PDOP:  From  King  et  at.  [1985,  Fig.  3.2|. 


Figure  4.  PDOP  versus  time 
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Hourly  observations  for  Monterey  were  obtained  from  the  NPS  Department  of 
Meteorology.  The  instruments  were  located  approximately  300  m  south  of  the 
Monterey  antenna  location  ( Figure  1 ). 


D.  PROCESSING  SOFTWARE 


A  complete  description  of  the  Trimble  supplied  Trimvec  software  can  be  found  in 
Trimble  Navigation  [l9S7b].  The  data  was  processed  using  the  Trimvec  Trim640  pro¬ 
gram.  Revision  AB.  T rim640  limits  processing  to  700  epochs,  so  the  first  700  epochs  for 
each  observing  session  are  used.  Sand  Point  was  used  as  the  reference  station  and  its 
coordinates  kept  fixed  in  the  least-squares  processing.  Sand  Point  was  chosen  as  the 
reference  station  because  four  satellite  availability  occurred  later  at  Sand  Point  than  at 
Monterey.  This  avoided  having  to  load  the  not-in-common  epochs  from  Monterey  at 
each  processing. 

Trim640  uses  the  C  A  code  derived  positions  obtained  at  the  lowest  PDOP  for  the 
initial  estimates  of  the  baseline  components.  Trim640  culls  the  best  C  A  code  position 
during  the  data  loading.  No  ionospheric  correction  is  provided  by  the  software,  and  only 
the  broadcast  ephemeris  can  be  used  to  compute  fixed  orbit  satellite  positions.  In  the 
triple  difference  processing,  the  only  parameters  estimated  by  the  least-squares  process¬ 
ing  are  the  baseline  components,  AX,  AY  and  AZ. 

A  modified  Hopfield  tropospheric  model  [Goad  and  Goodman,  1974]  is  used  to 
correct  the  carrier  phase  delay  caused  by  the  troposphere.  The  correction  is  a  function 
of  the  atmospheric  refractivity  computed  from  surface  meteorological  values  of  pressure, 
temperature,  and  humidity,  and  the  elevation  angle  of  the  satellite.  Larger  corrections 
are  required  for  low  elevation  angles,  as  the  signal  travels  a  longer  path  through  the 
troposphere.  The  model  corrects  for  at  least  90%  of  the  tropospheric  delay  [Remondi, 
198  h.  Tri/n640  allows  only  one  pressure,  temperature  and  humidity  entry'  for  each  site 


per  session. 


E.  PROCESSING  PROCEDURES 


1.  General 


To  study  optimized  times,  the  data  from  each  observing  session  were  segmented. 
Each  successive  segment  contained  10  minutes  more  data  than  its  predecessor.  For  ex¬ 
ample.  for  an  entire  observing  session  that  started  when  four  satellites  were  available  and 
stopped  when  less  than  four  were  available  provides  175  minutes  of  observations.  The 
first  segment  will  use  the  first  10  minutes  of  data  ,  the  sixth  segment  will  use  the  fust  60 
minutes,  while  the  eighteenth  will  use  all  175  nunutes.  For  each  segment,  the  entire 
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processing  u'as  restarted  from  the  data  loading.  Reloading  the  data  for  eacii  segment 
takes  considerably  longer  than  using  the  Trim640  option  to  Hag  data  for  processing,  hut 
reloading  was  done  so  that  the  processing  does  not  use  a  best  C  A  code  position  from 
later  in  the  observing  session. 

Convergence  of  the  least-squares  solution  was  achieved  by  doing  five  iterations 
using  every  tenth  triple  difference  formed  from  every  tenth  double  difference,  followed 
by  five  iterations  decreasing  the  triple  and  double  dilference  increments  to  five,  and 
finally  five  iterations  using  all  triple  differences  formed  from  all  double  differences. 

Trim640  rejects  those  observations  whose  residuals  exceed  a  multiple  of  the 
mean  residual.  The  multiple  of  the  mean  residual  is  known  as  the  edit  multiplier. 
Trim640  uses  3.5  as  the  default  value  for  the  edit  multiplier.  1  used  the  default  value  for 
the  initial  processing. 

Any  segment  that  had  more  than  ten  percent  of  its  observations  rejected  or 
whose  solution  slope  distance  standard  deviation  (a,)  was  greater  than  10  m  was  re¬ 
processed.  The  reprocessing  was  identical  to  the  initial  processing  except  that  before 
invoking  the  triple  difference  process  the  pseudoranges  for  both  stations  were  subjected 
to  separate  least-squares  adjustment.  The  pseudorange  processing  improves  the  C  A 
code  derived  initial  estimates  for  the  baseline  components  and  corrects  the  carrier  beat 
phase  time  tags.  The  carrier  beat  phase  time  tags  are  computed  from  the  C  A  code 
times,  and  are  earlier  than  the  C  A  code  times. 

If  the  pseudorange  processing  failed  to  lower  the  rejections  to  ten  percent,  the 
edit  multiplier  was  increased  until  the  rejections  reached  ten  percent.  A  ten  percent  re¬ 
jection  level  was  observed  for  a  few  sessions  and  always  occurred  within  the  first  thirty 
minutes  of  observations.  The  data  were  transferred  to  the  Naval  Postgraduate  School's 
IBM  3033  computer  for  analysis.  The  data  were  analysed  and  graphics  produced  using 
the  A  PL- based  GRAFSTAT  program. 

To  study  the  effects  of  reducing  observation  tune,  five  case  studies  were  under¬ 
taken  in  which  the  observation  start  times  were  changed  for  processing.  Each  case  study 
followed  the  processing  procedures  outlined  above. 

2.  Batch  Processing 

Processing  is  performed  in  a  batch  mode.  A  batch  file  passes  parameters  to  a 
template.  Trimble  supplies  command  files  that  tell  the  Trim640  to  use  the  template  pa¬ 
rameters  in  processing. 
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The  batch  file  is  built  using  the  program  Bath'd  (Appendix  At.  BaihUi  builds  a 
batch  file  by  providing  the  appropriate  tile  names  and  start  and  stop  times.  BaihlJ 
computes  the  appropriate  meteorological  parameters  lor  each  segment  by  locating  the 
applicable  weather  observations  from  the  weather  observation  tile,  interpolating  values 
at  the  start  time,  computing  running  means  from  each  hourly  weather  observation,  then 
interpolating  the  running  means  to  the  stop  time  of  each  segment.  Two  millibars  were 
subtracted  from  the  SLA-TAC  sea-level  pressure  to  compensate  for  the  20-m  elevation 
above  sea-level  for  the  Sand  Point  antenna. 

Initially,  processing  was  done  on  an  IBM  XT  with  a  math  coprocessor  and  a 
hard  disk.  Processing  the  18  segments  of  an  observing  session  took  ten  hours  of  com¬ 
puting  time.  Later,  processing  was  performed  on  an  80286  based  microcomputer  run¬ 
ning  at  10  MHz  with  an  S028'-8  math  coprocessor  that  reduced  the  processing  time  to 
three  hours. 

Two  minor  problems  with  Trim640  were  discovered  during  the  processing. 
First,  large  values  in  range  differences  were  found  when  using  the  range  differences 
rather  than  the  pseudoranges  to  improve  the  C  A  code  positions.  The  data  were  for¬ 
warded  to  Trimble  Navigation  for  evaluation  and  an  error  was  found  in  their  software. 

I  he  error  had  no  apparent  effect  upon  carrier  phase  difference  processing.  Second. 
Tnm64i)  is  incompatible  with  one  or  more  of  the  AST  Research.  Inc.  device  drivers 
supplied  with  the  80286  microcomputer:  ram  disk,  print  spooler,  and  extended  memory. 
Removing  the  drivers  allowed  Trim640  to  execute  normally. 

At  the  conclusion  of  the  batch  processing,  the  slope  distance,  the  baseline 
components,  their  standard  deviations  (a,),  the  number  of  observations,  the  number  of 
observations  rejected,  and  the  RMS  cycle  fits  were  extracted  from  the  Trim640  output 
file  and  collected  into  files  that  held  the  data  for  a  particular  segment  for  each  case  study. 

F.  ANALYSIS  PARAMETERS 

The  statistical  parameters  that  will  be  used  to  evaluate  the  results  are  the  error,  the 
sample  mean  and  the  standard  deviation  defined  as  [Davis  ei  al ..  1981): 


Et  Error  for  Jay  d 
C ,  Measured  component  for  day  d 

f'j  Expected  values 
E  Mean  error  or  bias 
a1  Variance 
a  Standard  deviation  of 

a i  Standard  deviation  of  E 

Accuracy  describes  the  closeness  between  the  measurements  and  the  expected  values 
(Davis  et  a/.,  1981).  The  degree  of  accuracy  is  determined  to  the  magnitude  of  the  mean 
error  (  E ).  The  repeatability  of  the  measurements  will  be  expressed  in  terms  of  2a  because 
it  approximates  the  95° o  confidence  level  for  single-dimension  measurements  (Federal 
Geodetic  Control  Committee,  1986|.  The  slope  distance  and  individual  baseline  com¬ 
ponents  are  one-dimensional  measurements. 

G.  DATA  AVAILABILITY 

Observations  were  made  simultaneously  at  Monterey  and  Sand  Point  for  an  eight- 
week  period  beginning  29  September  1987.  Observations  were  made  Tuesday  through 
Saturday  except  the  days  after  federal  holidays.  Forty  observing  sessions  were  con¬ 
ducted  of  which  2S  were  used  in  the  analysis  and  are  listed  in  Table  5.  The  remaining 
12  dnss  of  observations  were  not  used  in  the  analysis  for  various  reasons,  which  are 
listed  in  Table  6. 

For  brevity  the  observing  days  will  be  referred  to  by  their  Julian  day.  Times  in 
Table  5  are  given  in  Pacific  Standard  Time  (PST)  rather  than  Universal  Coordinated 
Time  (UTC)  for  ease  in  the  later  discussions  on  diurnal  effects. 
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Table  6.  OBSERVATION  DAYS  NOT  USED  IN  THE  ANALYSIS. 


Date 

Reason  for  not  analyzing 

10  15  S7 

No  SEA-TAC  weather  observations 

10  16  87 

No  satellites  at  Monterev  for  first  lo  minutes 

10  23  87 

Slope  distance  a,  >  10  m  for  first  10  minutes 

10  30  87 

Disk  error 

10  31  87 

No  SEA-TAC  weather  observations 

1 1  06  87 

Unhealthv  satellites 

11  07  87 

Unhealthv  satellites 

11  17  V 

Unhealthy  satellites 

11  18  87 

Unhealthy  satellites 

11  19  87 

Unhealthv  satellites 

1 1  20  87 

Unhealthv  satellites 

11  24  87 

No  satellites  at  Monterey  for  first  20  minutes 

Day  283  was  processed  with  Monterey  as  the  fixed  reference  station  because  the 
data  set  would  not  partition  into  lo-minute  segments  when  Sand  Point  was  used  as  the 
reference  station. 

The  da\s  that  were  missing  one  or  two  segments  were  excluded  from  the  analysis, 
so  that  changes  m  the  sample  variances  would  not  be  due  to  unequal  sample  populations 
between  the  segments. 

Reprocessing  the  first  segment  for  23  October  failed  to  reduce  the  slope  distance 
sigma  to  less  than  In  m  because  that  segment  had  only  nine  triple  difference  observa¬ 
tions  because  the  receiver  frequently  lost  lock  on  the  satellites.  That  day  was  not  used 
in  the  analysis,  even  though  its  other  segments  had  slope  distance  sigma's  less  than  10 
m  without  reprocessing. 
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IV.  RESULTS  AND  DISCLSSION 


A.  GENERAL 


To  study  optimized  times,  five  cases  are  studied: 

1  Process  all  data  when  four  satellites  at  least  15'  above  horizon 

2  Begin  processing  when  live  satellites  at  least  15°  above  the  horizon 

3  Begin  processing  as  m  2,  but  delete  fifth  satellite 

4  Begin  processing  data  40  minutes  later  than  in  1 

5  Begin  processing  data  70  minutes  later  than  in  1 

Case  1  is  essentially  the  processing  of  the  full  data  set.  Each  of  the  other  cases  is  a 
subset  of  Case  1.  Trim64()  allows  the  user  to  designate  at  what  time  within  the  full  ob¬ 
serving  period  that  the  data  loading  should  begin.  For  Case  3,  Trun640' s  ability  to  (lag 
data  was  used  to  exclude  the  fifth  satellite  (SV  12)  from  the  processing.  Each  day  of 
Table  5  was  processed  for  each  of  the  five  cases. 

During  the  course  of  the  discussions  it  will  be  necessary  to  distinguish  between  ob¬ 
servation  periods  and  the  time  of  the  observations  fixed  with  respect  to  satellite  geom¬ 
etry.  As  the  satellite  geometry  (PDOP)  begins  with  the  availability  of  four  satellites,  the 
time  of  observations  can  be  defined  in  terms  of  the  Case  1  start  time.  Observation  pe¬ 
riods  are  determined  from  the  start  time  for  each  case.  Times  of  observations  will  given 
as  equivalent  Case  1  times  and  is  obtained  by  adding  the  case  observation  length  to  the 
case  s  time  ofTset  from  Case  1  (Case  2  and  3's  offsets  are  20  minutes;  Case  4.  40  minutes 
and  Case  5.  "0  minutes). 

B.  ACCURACY 

1.  Slope  Distance 

The  slope  distance  errors  for  Cases  1.  2.  4.  and  5  are  presented  in  Table  7.  Then 
results  show  that  accuracy  to  better  than  1.0  ppm  is  achieved  for  any  observation  period, 
but  that  there  are  differences  among  the  cases  and  with  changing  observation  periods 
within  each  case. 

Cases  1,  2.  and  4  exhibit  similar  behavior  as  the  observation  period  increases  - 
they  become  less  positive  (or  more  negative)  as  the  observation  period  increases  until 
they  reach  a  minimum,  then  they  reverse  their  trends  and  become  less  negative.  Positive 
error  indicates  that  the  measured  baseline  is  longer  than  the  standard  values,  so  that  the 
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more  observations  are  included  in  the  solutions.  While  the  minimum 


occurs  alter  U, 


ferent  observing  periods  <  14t>.  110.  mo  minutes  for  Cases  1,  2.  and  4  respectively  i.  thev 
occur  at  about  the  same  absolute  time  with  respect  to  the  Case  1  start  time  <  I4n.  1’". 
and  130  minutes  for  Cases  1.2.  and  4  respectively).  As  the  start  times  occur  later  with 
each  case,  the  errors  for  shorter  observing  periods  become  less  positive  and  the  ranges 
of  the  errors  lor  each  case  decrease  (the  range  of  Case  1  errors  is  3o  to  -IX. 4  cm;  -  ~.S 
to  -IX. 1  cm.  lor  Case  2;  and  -10.3  to  -15.1  cm.  lor  Case  4).  The  error  after  the  entire 
observation  session  decreases  with  later  observing  start  times  t  -14. S  cm.  -V  I  cm.  and 
-0.4  cm  for  Cases  1.  2.  and  4  respectively). 

Case  5  behaves  similarly  to  the  previous  cases  except  that  its  minimum  occurs 
after  only  10  minutes  of  observations  and  adding  more  observations  causes  the  error  to 
become  positive.  The  error  is  largest  after  the  entire  observing  session  ( 3S.4  cm).  Cases 
1.  2,  and  4  start  their  observations  prior  to  the  first  PDOP  peak  that  occurs  at  60  min¬ 
utes  (figure  4)  while  Case  5  starts  after  the  PDOP  peak,  so  that  starting  the  observa¬ 
tions  close  to  the  larger  PDOP  peak  and  including  the  PDOP  peak  observations  can 
reduce  the  error  and  the  required  observation  period.  The  effect  of  the  PDOP  peak  upon 
the  mean  slope  distance  error  is  readily  apparent  in  comparing  Cases  4  and  5.  Case  4 
remains  negative  without  the  early  positive  error,  and  Case  5  remains  positive  as  it  lacks 
most  of  the  observations  from  about  the  first  PDOP  peak. 

fhe  first  PDOP  peak  at  60  minutes  differs  from  the  second  PDOP  peak  in  that 
it  has  a  higher  value,  is  symmetric,  and  occurs  farther  from  the  four  satellite  observation 
start  time  than  the  second  peak  is  from  the  four  satellite  observation  stop  time,  the 
second  peak  occurs  with  lower  elevation  angles  than  the  first  which  implies  larger 
tropospheric  and  ionospheric  errors. 

Ihe  results  for  Case  3  (excluding  the  fifth  satellite)  slope  distance  errors  are 
prevented  in  Table  8.  Case  3  was  studied  because  the  Case  2  results  showed  a  decrease 
in  the  initial  slope  distance  errors,  and  Case  3  was  to  study  the  effects  of  observing  five 
satellites.  The  results  show  that  using  only  four  satellites  makes  a  difference  of  only  a 
few  centimeters  from  the  results  using  all  available  satellites.  It  should  be  noted  that 
Case  3  uses  only  three  satellites  once  SV  6  sets  after  100  minutes  of  observations.  The 
Case  3  minimum,  -11.5  cm,  was  less  neeative  than  the  Case  2  minimum  of  -18.1  cm. 
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2.  Baseline  Components 

1  he  AX  and  AY  accuracies  are  better  than  lo.it  ppm  for  all  obvervine  periods 
while  AZ  accuracy  is  better  tiian  1.0  ppm  lor  ail  observing  periods.  The  accuracv  of  the 
baseline  components  is  expected  to  be  less  than  the  accuracy  of  the  baseline  because  the 
baseline  errors  are  mostly  perpendicular  to  the  baseline  itself  [Remondi.  19s4j.  The 
baseline  component  results  (Tables  9,  |o.  and  11)  show  that  the  AX  and  AY  errors  are 
greater  than  the  baseline  components  while  the  AZ  errors  are  about  the  same  order  of 
magnitude  as  the  slope  distance  enors.  The  AX  and  Ah'  errors  are  negatively  correlated 
which  is  the  result  of  the  correlations  of  the  triple  dilTerences  and  the  semi-circular  tracks 
of  the  satellites  (Figure  5).  Case  3  (Table  Si  shows  little  difference  from  Case  2  in  the 
AX  error,  and  a  more  negative  Ah'  error  is  olTset  by  a  less  negative  AZ  error. 

The  smallest  mean  errors  for  the  baseline  components  are  found  at  various  ob¬ 
serving  periods.  Zero  mean  error  for  all  the  baseline  components  is  achieved  with  fewer 
observations  as  the  observing  start  times  occur  closer  to  and  before  the  larger  PDOP 
peak. 

Table  9.  MEAN  AX  ERROR  FOR  CASES  1.  2.  4.  AND  5 


Observation 

Case  1 

Case  2 

-  - 

Case  4 

Cave  5 

Period 

[mini 

cm 

ppm 

cm 

ppm 

cm 

ppm 

cm 

ppm 

10 

-252.4 

-6. 1 3 

-296.3 

-7.19 

-136.3 

-3.30 

29.6 

0.73 

2o 

-255.4 

-6.21 

-256.0 

-6.22 

-136.4 

-3.30 

83.9 

2.04 

3o 

-255.5 

-6.21 

-212.1 

-5.15 

-107.6 

-2.62 

132.9 

3.23 

do 

-238.0 

-5.78 

-194.8 

-4.74 

-80. 1 

-1.94 

183.9 

4.47 

5o 

-211.2 

-5.12 

- 1 6  / .  8 

-4.08 

-52.7 

-1.29 

191.7 

4. no 

00 

-197.0 

-4.79 

-137.4 

-3.33 

-24.8 

-0.61 

192.9 

4.69 

"o 

-171. S 

-4.17 

-107.3 

-2.60 

21.9 

0.53 

193.4 

4.70 

so 

-145.1 

-3.52 

-76.8 

-1.S7 

48.4 

1.17 

195.4 

4.75 

90 

-116.9 

-2.84 

-31.8 

-0.7S 

67.7 

1 .65 

201.7 

4.91 

100 

-87.3 

-2.11 

3.1 

0.08 

82.3 

1.99 

208.7 

5.08 

1  10 

-45.2 

-1.09 

28.4 

0.70 

96.2 

2.33 

221.0 

5.37 

120 

-10.7 

-0.27 

48.6 

1.18 

112.0 

2.72 

- 

130 

15.2 

0.36 

65.7 

1.60 

124.8 

3.04 

- 

140 

35.6 

0.87 

83.3 

2.02 

139.9 

3.40 

. 

150 

53.6 

1.30 

97.7 

2.37 

- 

- 

- 

160 

70.7 

1.72 

110.1 

2.67 

- 

- 

170 

84.3 

2.04 

. 

. 

. 

• 

_ 

175 

91.7 

2.24 

- 

" 

- 

- 

Table  10.  MEAN  AY  ERROR  FOR  CASES  I.  2.  4.  AND  5 


Observation 

Period 

|min| 

Case  1 

Case  2 

Case  4 

Case  5 

cm 

ppm 

cm 

ppm 

cm 

lo 

1SS.0 

2.63 

199.2 

2.“S 

1 10.0 

1.54 

9.1 

0.13 

2'  > 

1S1.2 

2.53 

181.4 

2.53 

106.3 

1.49 

-28.0 

-11.39 

30 

IT9.3 

2.50 

152.4 

2. 1  3 

97.6 

1.31 

-68.6 

-0  96 

40 

171.2 

2.39 

139.3 

1.95 

79.  s 

111 

-95.2 

-1.33 

50 

153.0 

2.14 

125.3 

1.75 

66.6 

0.93 

-92.6 

-1.29 

60 

141.0 

1.97 

109.5 

1.53 

52.2 

o.~3 

-91.9 

- 1 .28 

“0 

127.7 

1.78 

94.3 

1.32 

29.3 

0.41 

-89." 

-1.25 

so 

1  14.0 

1.59 

79.3 

1.11 

21.6 

0.30 

-89. 3 

-1.25 

00 

100.5 

1 .40 

5S.4 

o.S2 

11.1 

0.16 

-95.6 

-1.33 

I'M) 

86.4 

1.21 

46.4 

0.65 

4.0 

0.06 

-105.4 

-1.4" 

110 

66.7 

0.93 

34.0 

0.47 

O.  J 

-0.o5 

-121.0 

-1.69 

120 

55  S 

0.78 

23  5 

0.33 

-15.S 

-0.22 

- 

- 

130 

43.3 

0.60 

14.8 

0.21 

-27.7 

-0.39 

- 

- 

140 

33.9 

0.47 

2.9 

o.o4 

-42.5 

-0.59 

- 

- 

15o 

25.5 

0.36 

-9.0 

-0.13 

- 

- 

- 

- 

160 

14.4 

0.20 

-17.0 

-0.24 

- 

- 

- 

- 

170 

4.0 

0.06 

- 

- 

- 

- 

- 

- 

1 T5 

-l.S 

-0.03 

- 

- 

- 

- 

- 

- 

Table  11.  MEAN  AZ  ERROR  FOR  CASES  I.  2.  4.  AND  5 


Observation 

Period 

[min| 

Case  1 

Case  2 

Case  4 

Case  5 

cm 

ppm 

cm 

ppm 

cm 

ppm 

10 

-S2.2 

-0.90 

-33.2 

-03.6 

-1 1.0 

-o  12 

-1".4 

-n  19 

20 

-68.0 

-0.75 

-30.2 

-0.33 

“7  *7 

- ' .  J 

-0 .08 

-23.6 

1.26 

30 

-54.8 

-0.60 

-23.4 

-0.36 

-8.5 

-0.1)9 

-'•2.0 

-li  'h 

40 

-47.3 

-0.52 

-17.0 

-0.19 

-8.9 

-0.  lo 

-34.  S 

-d.3S 

50 

-38.1 

-0.42 

-14.6 

-0. 1 6 

-11.1 

-0. 1 2 

-o  41 

60 

-28.9 

-0.32 

-12.3 

-0.13 

-15.1 

-0.17 

-39.1 

-0.42 

7o 

-23.8 

-0.26 

-11.6 

-0.13 

-15.1 

-0.17 

-44." 

-i  i  4a 

SO 

-20.0 

-0.22 

-13.1 

-0.14 

-18.5 

-0.20 

-49  5 

a)  >4 

90 

-18.1 

-0.20 

-12.4 

-0.14 

-18.9 

-0.20 

.<2.7 

-0.58 

100 

-17.9 

-0.20 

-15.1 

-0.17 

*.  «-  •  4m 

-o.24 

-54.2 

.1  1. 

110 

-15.6 

-0.17 

-15.1 

-0.17 

-25.7 

-0.2S 

-56.5 

-i  i.o2 

120 

-17.5 

-0.19 

-17.2 

-0.19 

-27.3 

-0.30 

- 

- 

130 

-16.4 

-0.18 

-19.S 

-0.22 

-28.0 

-0.31 

- 

- 

140 

-17.9 

-0.20 

-21.0 

-0.23 

-29.3 

-0.32 

- 

- 

150 

-19.9 

-0.22 

-21.3 

-0.23 

- 

- 

- 

- 

160 

-20.3 

-0.22 

-21.9 

-0.24 

- 

- 

- 

- 

170 

-20.1 

-0.22 

- 

- 

- 

- 

- 

- 

175 

-20.1 

-0.22 

- 

- 

- 

- 

- 

- 
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C.  REPEATABILITY 
I.  Slope  Distance 

The  day-to-day  repeatabilities,  represented  as  the  2a  level,  lor  Cases  1,  2,  4,  and 
5  are  presented  in  Table  12,  and  for  Case  3,  in  Iable  13.  All  the  cases  achieve  l.o  ppm 
repeatability  for  any  observing  period  except  Case  1  which  requires  20  minutes  of  ob¬ 
servations.  Repeatability  eventually  reaches  better  than  0.5  ppm  after  60  minutes  of 
observations  lor  any  case. 

The  minimum  2 a  levels  for  all  cases  are  reached  at  the  SO  to  90  niinute  time  of 
observation,  which  is  about  30  minutes  after  the  larger  PDOP  peak.  A  slight  increase 
in  the  2 a  level  is  centered  about  the  120  to  130  minute  time  of  observation  for  all  cases 
which  is  near  the  PDOP  minimum,  after  which  the  2a  level  decreases  slightly  as  obser¬ 
vations  from  the  second  PDOP  peak  are  included  in  the  solutions.  Case  4  had  the 
narrowest  range  of  2cr  values,  and  Case  1  had  the  widest  range  of  2cr  values. 


Table  12.  SLOPE  DISTANCE  2-SIGMA  VALUES  FOR  CASES  1.  2.  4,  AND  5 


Observation 

Period 

(min] 

Case  1 

Case  2 

Case  4 

Case  5 

cm 

ppm 

cm 

ppm 

cm 

ppm 

cm 

ppm 

■HI  - 

1.44 

64.8 

0.53 

5^.8 

52.0 

0.42 

0.96 

63.2 

0.51 

56.8 

'  s  ■ 

53.4 

0.43 

84.1 

0.68 

58.6 

0.4S 

59.2 

1  1  1 

60.8 

0.49 

.'VH  ' 

69.8 

0.57 

58.4 

0.47 

53.8 

57.0 

0.46 

50 

64.0 

0.52 

58.6 

0.48 

54.4 

0.44 

55.4 

0.45 

60 

61.8 

0.50 

55.2 

0.45 

57.6 

0.47 

55.0 

0.45 

70 

60.4 

0.49 

55.0 

0.44 

59.6 

0.48 

53.6 

0.44 

80 

57.4 

0.47 

56.8 

0.46 

60.4 

0.49 

53.6 

0.44 

90 

56.2 

0.46 

59.0 

0.48 

60.8 

0.49 

57.6 

0.46 

100 

57.2 

0.47 

60.0 

0.49 

59.6 

0.48 

62.6 

0.51 

no 

59.2 

0.48 

60.4 

0.49 

57.S 

0.47 

67.8 

0.55 

120 

60.0 

0.49 

59.6 

0.48 

56.2 

0.46 

130 

61.0 

0.50 

59.0 

0.48 

56.0 

0.46 

140 

60.6 

0.49 

58.2 

0.47 

55.4 

0.45 

150 

60.0 

0.49 

57.8 

- 

- 

- 

- 

59.0 

0.48 

58.4 

- 

. 

. 

58.6 

0.48 

- 

- 

. 

58.6 

0.48 

- 

- 

- 
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Table  13.  CASE  3:  2-SIGMA  VALl'ES 


Observation 

Period 

Imin] 

Slope 

Distance 

AX 

AY 

AZ 

cm 

ppm 

cm 

ppm 

cm 

ppm 

cm 

ppm 

H» 

62.4 

0.51 

326.0 

7.92 

140.8 

1.97 

7  6 . 6 

O.S4 

2') 

59.4 

0.49 

304.2 

7.3  6 

135.8 

1.90 

80.  S 

( 1.89 

30 

57.8 

0.47 

318.8 

7.75 

169.0 

2.36 

76.6 

o.84 

40 

58.0 

0.47 

3ol.O 

7.31 

1 73.4 

2.42 

82.2 

0.90 

50 

58.8 

0.48 

287.4 

6.98 

195. S 

2.60 

95.4 

1.04 

54.0 

0.44 

->7->  7 

6.61 

197.0 

2.75 

97.0 

1.06 

70 

53.0 

0.43 

251.8 

6.12 

199.8 

2.79 

95.8 

1.0  5 

SO 

53.6 

0.44 

244.6 

5.94 

199.4 

2.79 

S9.0 

0.98 

00 

56.4 

0.46 

235.6 

5.72 

187.0 

2.61 

83.4 

0.91 

100 

58.0 

0.47 

219.8 

5.34 

165.8 

2.32 

76.2 

0.S4 

1 10 

59.2 

0.48 

199.2 

4.84 

13S.8 

1.94 

71.8 

o,79 

120 

59.0 

0.4S 

188.8 

4.59 

127.4 

1.78 

6S.2 

0.75 

130 

59.0 

0.48 

184.2 

4.48 

122.6 

1.71 

65.0 

0.71 

140 

58.8 

0.48 

ISS.4 

4.5S 

127.2 

1.78 

65.4 

0.72 

150 

58. S 

0.48 

190.4 

4.62 

128.4 

1.79 

64.8 

0.71 

160 

5".  8 

0.47 

20  3.0 

4.93 

137.0 

1.91 

63.8 

0.70 

2.  Baseline  Components 

All  the  baseline  components  have  repeatabilities  better  than  10.0  ppm  for  any 
observing  period  except  for  the  Case  1  AX,  which  required  20  minutes  of  observations 
(Tables  13,  14.  15.  and  16). 

It  is  interesting  to  note  that  while  the  AX  2<x  values  for  the  first  segment  of 
Cases  2,  3,  4,  and  5  are  less  than  Case  I  s  first  segment,  the  final  Case  1  2«r  value  is  less 
than  the  final  segment  of  any  other  Case.  The  final  2cr  values  for  Cases  2,  3,  4,  and  5 
are  greater  than  the  Case  1  2<7  values  after  an  equivalent  number  of  observations. 

The  minimum  2<7  levels  for  AX  and  AY  occur  at  or  near  the  end  of  the  observing 
sessions  for  Cases  1,2,  4  and  5.  For  Case  3  the  minimum  2 a  levels  occur  after  130  min¬ 
utes  of  observations.  The  minimum  A Z  2<j  level  occurs  after  various  observation  peri¬ 
ods,  but  generally  in  the  vicinity  of  130  to  150  minute  observation  time,  which  is  between 
the  PDOP  minimum  and  the  second  PDOP  peak.  Case  3  behaves  in  an  opposite  fashion 
from  the  other  cases  in  that  its  minimum  AX  and  AY  2 a  levels  occur  at  the  150-minute 
time  of  observation  while  its  A Z  minimum  occurs  at  the  end  of  the  observation  period. 
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Table  14.  AX  2-SIGMA  VALUES  FOR  CASES  1.  2,  4.  AND  5 


Observation 

Period 

[min] 

Case  1 

Case  2 

Case  4 

Ca'-e  5 

cm 

ppm 

cm 

ppm 

ppm 

cm 

ppm 

10 

646.4 

15.69 

3JC2 

8.45 

422.0 

10.27 

39  V  2 

9.1.0 

:<) 

531.4 

12.93 

311.4 

7.5" 

343.4 

8.34 

3  "5.4 

9.12 

30 

345.0 

8.41 

312.2 

".58 

312  2 

"32 

3  "5.0 

8.8" 

40 

248.2 

6.03 

3 1 1 .0 

7.31 

286. 8 

6.97 

33<v2 

8.1" 

5o 

239.6 

5.83 

286.0 

6.95 

251.0 

6.10 

315.8 

".43 

60 

230.4 

.5.59 

209.2 

6.54 

24o.S 

6.00 

2"2.S 

(.  (>3 

"0 

221.2 

5.39 

246.6 

5.99 

233.4 

5.67 

252.8 

(v  14 

So 

209.2 

5.10 

243.2 

5.91 

212.8 

5.1" 

231.0 

5  ol 

90 

195.0 

4.74 

221.4 

5.38 

188.4 

4.43 

224.4 

5.45 

100 

202.0 

4.91 

194.8 

4.73 

168.8 

4.10 

206  6 

5.o2 

110 

180.8 

4.39 

168.8 

4.10 

159.2 

3.S7 

203.4 

4.94 

120 

162.4 

3.94 

155.2 

3.77 

155.6 

3. "8 

- 

130 

141.4 

3.44 

146.4 

3.56 

141.4 

3.43 

- 

140 

134.6 

3.26 

144.4 

3.51 

139.8 

3.40 

- 

150 

129.0 

3.11 

133.0 

3.23 

- 

- 

- 

160 

129.4 

3.14 

134.4 

3.27 

- 

- 

- 

170 

121.8 

2.96 

- 

- 

- 

- 

- 

175 

121. S 

2.96 

- 

- 

- 

- 

- 

Table  15.  AY  2-SIGMA  VALUES  FOR  CASES  I.  2.  4.  AND  5 


Observation 

Period 

[min| 

Case  1 

Case  2 

Case  4 

Case  5 

cm 

ppm 

cm 

ppm 

cm 

ppm 

cm 

ppm 

10 

110.8 

1.55 

204.0 

2.85 

244.2 

3.41 

325.0 

4.54 

20 

130.4 

1.82 

2.67 

224.2 

3.13 

311.4 

4.34 

30 

140.8 

1.97 

190.6 

2.67 

3.08 

29S.4 

4.17 

130.8 

1.83 

193.2 

2.70 

221.0 

3.09 

264.8 

3.70 

135.4 

1.89 

192.2 

2.69 

212.0 

2.96 

237.0 

3.31 

142.8 

1.99 

188.6 

2.63 

3.00 

200.0 

2.79 

70 

144.4 

2.02 

185.2 

2.59 

214.6 

2.S5 

181.4 

2.53 

80 

144.6 

2.02 

185.0 

2.58 

181.6 

2.54 

170.2 

2.38 

90 

145.4 

2.03 

173.4 

2.42 

147.4 

2.06 

177.2 

2.4S 

too 

145.0 

2.10 

152.4 

2.13 

133.8 

1.87 

174.6 

2.44 

110 

139.0 

1.94 

127.6 

1.7S 

124.6 

1.74 

182.2 

2.55 

120 

126.8 

1.77 

115.2 

1.61 

124.8 

1.74 

- 

130 

108.0 

1.51 

110.6 

1.55 

119.0 

1.66 

- 

140 

102.6 

1.43 

112.6 

1.57 

120.4 

1.68 

- 

150 

100.0 

1.40 

109.0 

1.52 

- 

- 

- 

160 

101.8 

1.42 

118.0 

1.65 

- 

- 

- 

170 

100.0 

1.40 

- 

- 

- 

- 

- 

175 

101.8 

1.42 

- 

• 

- 

- 

- 

Table  16.  AZ  2-SIGMA  VALUES  FOR  CASES  I.  2,  4.  AND  5 


Case  1 


ppm 

cm 

PPm 

cm 

■dr®. 

65.2 

75.4 

0.76 

0.81 

92. S 

8'. 4 

(1.54 

69.  S 

0.7' 

S2.S 

0.60 

71.6 

u.~9 

78.4 

0.61 

71.6 

0.79 

72.o 

0.66 

68.4 

o.  -r5 

68.4 

0.69 

63.8 

0.7() 

65.0 

0.07 

59.  S 

0.70 

59. S 

0.65 

58.2 

0.64 

55.2 

0.64 

55.4 

0.6 1 

53.2 

0.62 

53.4 

0.59 

52.2 

0.61 

52.8 

0.5S 

53.6 

0.60 

52.4 

0.5S 

54.0 

0.59 

54.2 

0.59 

53.6 

0.59 

54.4 

0.60 

. 

0.61 

53.4 

0.59 

- 

0.62 

- 

- 

- 

0.61 

- 

- 

- 
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3.  Standard  Deviation  of  the  Mean 


The  repeatability  values  can  be  used  to  estimate  the  standard  deviations  of  the 
mean  errors  given  in  Tables  7  through  1 1  by  using  Equation  (5.2).  The  values  of  Tables 
12  through  16  should  be  divided  by  N  2S  (where  28  is  the  sample  population)  to  compute 
the  standard  deviations  of  the  means  (at  the  2 a  level).  Generally,  the  repeatabilities  w  ere 
about  five  times  the  magnitudes  of  the  mean  errors;  therefore,  the  uncertainties  of  the 
mean  errors  are  on  the  order  of  the  mean  errors  themselves.  Allowing  for  the  0.1  ppm 
uncertainty  in  the  baseline  and  the  baseline  components  and  for  the  possible  standard 
deviation  of  the  means,  accuracies  to  better  than  1.0  ppm  for  the  slope  distances  and 
10.0  ppm  for  the  baseline  components  remain  valid. 

D.  ERROR  EFFECTS 
1.  7-Day  Means 

Because  of  the  observations  were  made  over  a  long  period  of  time.  Case  1  was 
subdivided  into  four  groups  comprised  of  seven  consecutive  observation  days  to  study 
trends  in  the  slope  distance  error  to  identify  the  contribution  of  various  error  sources  to 
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optimizing  obsersing  times.  The  results  for  the  slope  distance  error  are  presented  in 


Table  17. 


Table  17.  CASE  I  ERROR:  7-DAY  MEANS 


The  seven-day  mean  slope  distance  errors  remain  below  1.0  ppm  for  all  groups 
for  any  observation  period.  The  differences  between  the  groups  are  mainly  in  the 
predominance  of  negative  errors  in  Groups  2  and  4  while  Groups  1  and  4  have  pre¬ 
dominantly  positive  errors. 

To  account  for  the  change  in  the  characters  of  the  groups,  the  change  in  the 
major  sources  of  error  are  examined. 

2.  Ephemeris  Errors 

Examining  the  AODE  (Age  of  Data  Ephemeris)  for  all  satellites  during  the 
eight-week  observing  period,  showed  that  the  oldest  AODE  for  a  single  satellite  was  nine 
hours,  but  most  were  less  than  five  hours.  Averaging  the  AODE  for  all  satellites  for 
each  observing  day  showed  a  range  of  seven  hours  during  the  eight  weeks  of  observa¬ 
tions  (Figure  6).  Assuming  12  hours  is  the  largest  AODE.  the  maximum  ephemeris 
linearization  error  would  be  50  cm.  Using  Equation  (2.3)  results  in  0.025  ppm  baseline 
relative  error. 
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Figure  6.  Mean  age  of  data  (ephemeris) 

The  changes  in  the  seven-day  mean  AODE's  do  not  correspond  well  to  the 
changes  in  the  groupings,  especially  as  the  largest  change  in  AODE  from  group  3  to  4 
does  not  correspond  to  a  similar  change  in  the  group  3  to  4  mean  slope  distance  errois. 
This  is  not  surprising  because  while  the  change  is  relatively  large,  the  magnitude  of  the 
orbital  errors  is  expected  to  be  small. 

The  ephemeris  error  will  appear  as  a  bias  during  any  one  observing  session.  By 
averaging  over  eig  .t  weeks  some  of  the  errors  will  cancel  and  appear  as  variance.  Be¬ 
cause  of  the  short  observing  sessions,  little  variation  in  the  ephemeris  error  is  expected 
during  any  single  sessions 

3.  Ionospheric  Errors 

All  the  observing  sessions  completed  well  before  1600  local  standard  time  and 
approached  the  normally  ionospheric  activity  minimum  of  0600  (Figure  7).  The  near 
north-south  orientation  of  the  Monterey-Sand  Point  baseline  places  both  stations  in  the 
same  time  zone,  so  that  ionospheric  errors  between  the  two  stations  will  be  correlated 
to  some  degree  and  reduced  in  the  differencing. 
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The  weekly  Preliminary  Report  and  Forecast  of  S<>iar  Geophyaeal  Ihua  |Spnce 
Environment  Services  Center.  I^S'J  for  'he  eight  weeks  of  the  observ aliens  described  the 
solar  and  geomagnetic  activity  as  generailv  quiet  or  low. 

Though  the  observing  period  Occurred  during  a  general  lull  in  ionospheric  ac¬ 
tivity.  there  are  two  ionospheric  elVects  that  can  be  eas,l\  examined:  the  diurnal  effect 
as  the  observing  period'  occur  closer  to  the  diurnal  ionospheric  actr. nv  nunuiuun  and 
'poiudic  effects  horn  electron  f.uxes.  I  lie  change  m  lonosplieric  range  error  for  tl;e 
eight-week  observing  period  was  computed  It  am  Henson  and  Collier  |19S<\  Equation 
i.'ij  using  t!ie  dillerence  in  Total  Electron  Content  ■  T 1  <  i  determined  from  Spilker  |  I'fSn. 
I  igure  1-1 1|  and  Henson  and  Collier  |  1486.  I  iguie  2|  Tor  the  mean  observation  tunes  for 
the  first  and  lust  groups.  The  maximum  change  in  ranee  error,  because  of  diurnal 
ionospheric  changes  from  the  first  to  the  last  week,  is  approximately  1.5  m  which  by 
Equation  (2.3i  is  equivalent  to  a  change  in  the  relative  baseline  error  of  only  <>.u"5  ppm. 

While  the  nil's  ppin  trend  is  not  discernible  because  of  the  0.1  ppm  uncertainty 
in  the  baseline  and  the  larger  uncertainties  in  the  measurements,  the  mean  number  of 
observations  i Table  5»  increased  with  each  later  group.  The  mean  number  of  triple 
difference  observations  available  were  1519  ±45,  1522  r  22.  ]ool  ^124.  and  I "26  ±1"4 
tor  the  first  through  the  last  groups  respectiv  ely.  The  increased  number  of  observations 
could  be  attributed  to  stronger  signai-to-noise  ratio  with  the  reduced  ionospheric 
dispersion  as  the  observations  approach  the  diurnal  ionospheric  activity  minimum. 

The  group  to  group  changes  did  not  show  an  apparent  trend,  but  alternated 
between  predominantly  negative  and  positive  mean  slope  distance  errors  indicating  a 
sporadic  ellect.  or  complex  interaction  between  the  error  sources.  To  examine  sporadic 
ionospheric  effects,  daily  Electron  Fluence  (EF)  from  the  weekly  Preliminary  Report  and 
Forecast  of  Solar  Geophysical  Data  (Space  Environmental  Services  Center,  198"],  and 
seven-day  means  of  the  EF  were  computed  (Figure  S).  Electron  fluence 
(  electrons  —  cm  2 -day-' —sr'1  )  is  the  daily  average  of  electron  flux  with  energies 
greater  than  2  Mev  as  measured  by  the  GOES-7  satellite.  High  values  of  the  seven-day 
mean  EF  corresponded  to  predominantly  positive  seven-day  mean  slope  distance  errors 
(Groups  1  and  3).  and  low  mean  EF  corresponded  to  negative  mean  slope  distance  error. 

4.  Tropospheric  Errors 

Refractivity.  ,V.  was  computed  using  the  meteorological  parameters  from  the 
time  closest  to  the  middle  of  each  observation  session  and  Remondi  [1984.  equations 
2.30  and  2.31].  Possible  indicators  of  the  influence  of  the  troposphere  upon  the  baseline 
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solution  are  the  diiTerer.ee  in  reiractn.it>.  A.V.  between  the  baseline  stations  i  I  igure  *' > 
and  the  mean  .V  of  both  stations  if  igure  I'm  The  Jtttcrcnce  car.  indicate  tire  lose!  .  i 
correlation  of  the  tropospheric  errors.  Correlated  tropospheric  errcts  are  reduced 
differencing.  The  mean  .V(.V)  will  indicate  the  magnitude  oi  the  reii.ia.cn  cl  me  carrier 
signal  to  'noth  sites. 

ii.e  changes  in  the  dillerence  in  the  se\en-da\  mean  .V  ( |  igure  d.d  ic  have 
an  apparent  correspondence  to  the  changes  in  the  seven-day  group  mean  slope  distance 
errors.  Instead  A  shows  a  pattern  similar  to  the  group  slope  distance  errors 
i  figure  10).  I Iigh  mean  refractivaty  corresponded  to  positive  mean  slope  distame  errors 
(Groups  I  and  3 1,  and  low  mean  refractivuv  corresponded  with  negative  mean  slope 
distance  errors.  Ihe  ditTerence  between  t!ie  high  and  c he  low  mean  A  are  small  and  do 
not  readily  account  lor  the  magnitude  ot  the  changes. 


Figure  10.  Mean  Monterey-Sand  Point  refractivity 
5.  7-Day  Repeatability 

The  Group  repeatabilities  could  not  be  given  a  general  characterization  as  the 
Group  means  were  because  the  was  no  consistent  pattern  to  the  Group-to-Group 
changes  in  their  repeatabilities.  For  completeness,  the  2cr  values  for  the  various  error 
sources  are  presented  in  Table  18,  and  the  Group  2<7  values  are  presented  in  Table  19. 
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Table  18.  GROUP  2-SIGMA  VALUES  FOR  THE  ERROR  SOURCES 


Error  Source 

Group  1 

Group  2 

Group  3 

Group  4 

AODE 

[hours] 

1.6 

2.0 

1.6 

2.4 

Observation  End  Time 
-  0600  PST 
[minutes] 

24.8 

45.0 

37.6 

56.8 

Electron  Fluence 
[electrons  —  enr 1  —day-'  —  sr~'] 

610.0 

59.6 

294.8 

21.4 

AN 

27.4 

21.2 

27.8 

23.8 

V 

15.2 

10.6 

9.0 

17.8 

Table  19.  CASE  I:  2-SIGMA  VALl'ES.  7-DAY  GROUPINGS 


Observation 
Period 
[  mini 

Group  1 

Group  2 

Group  3 

Group  4 

cm 

ppm 

cm 

ppm 

cm 

ppm 

in 

1 20. S 

1 .06 

8 12 

— 

0.66 

155.4 

1.26 

— 

:;o.o 

— 

1  .so 

2o 

SO. 4 

0.73 

1  lo  2 

0.00 

141.8 

1.15 

124  S 

1  01 

.’o 

61.8 

0.50 

64.0 

o.52 

loS.2 

0  S8 

80.6 

0  ~3 

40 

t'O.S 

0.40 

5o.2 

o.46 

SO.  6 

O.'o 

64. S 

o53 

50 

58.0 

0.47 

53.6 

0.44 

'2.6 

0.59 

52.4 

0.43 

no 

59.8 

0  40 

53.4 

0.43 

6  V4 

0.53 

45  6 

0  3~ 

To 

62.0 

0.50 

46.2 

o.3S 

60.2 

0.49 

44.8 

0.36 

So 

57.0 

0.46 

40.0 

0.40 

53.4 

0.43 

4S.2 

0.39 

00 

5  AO 

0.46 

47.4 

0.30 

51.4 

0.42 

4'. 6 

0.39 

lot) 

55.6 

0.45 

4o.O 

0.37 

5o.4 

0.41 

50.0 

0.41 

1  10 

57. 4 

0.47 

46.8 

0.38 

5o.4 

0.41 

53.2 

0.43 

120 

61.8 

0.50 

46.8 

0.38 

48.2 

0.39 

55.4 

0.45 

13o 

64.6 

0.53 

44.4 

0.36 

45.4 

0.37 

58.0 

0.47 

140 

6  .5  0 

0.53 

43.0 

0.35 

43.4 

0.35 

60.2 

0,49 

15o 

64.4 

0.53 

43.4 

0.35 

42.2 

0.34 

60.2 

0.49 

160 

63.4 

0.52 

42. S 

0.35 

42.6 

0.35 

59.2 

0.48 

ro 

62.8 

0.51 

42.6 

0.35 

43.8 

0.36 

58.6 

0.48 

175 

63.0 

0  51 

41.6 

0.34 

45.0 

0.37 

57.6 

0.47 

E.  EFFECTS  OF  THE  C/A  CODE 

The  best  C  A  code  position  for  Monterey  was  extracted  from  each  observing  session 
for  comparison  with  the  carrier  phase  results.  The  coordinates  of  Sand  Point  were  sub¬ 
tracted  from  the  28-day  mean  of  the  Monterey  C  A  code  coordinates  to  form  the 
baseline  components.  The  slope  distances  were  computed  from  the  square  root  of  the 
sum  of  the  squares  of  the  baseline  components.  The  mean  errors  of  the  C  A  code  de¬ 
rived  baseline  were  on  the  order  of  the  carrier  phase  mean  errors,  but  the  2 a  levels  are 
several  times  greater  than  the  carrier  phase  results  (Table  20).  These  results  are  for  one 
instantaneous  position  determination  and  would  improve  with  time  averaging  within  the 
observing  period  or  Doppler-smoothing. 

The  C  A  code  has  its  most  significant  impact  upon  the  differenced  carrier  phase 
solutions  when  there  are  few'  carrier  phase  observations.  Table  21  shows  that  reproc¬ 
essing  the  carrier  phase  observations  was  required  for  several  segments  because  either 
more  than  ten  percent  of  the  observations  were  rejected  or  the  slope  distance  solution 
sigma  was  greater  than  10  m,  and  that  all  the  reprocessed  segments  were  confined  to  the 
first  two  segments. 


Table  20.  BEST  C/A  CODE  RCSl'LTS  AND  ERRORS 


Component 

Distance  •  m> 

Reiatise  1. 1  tor  i  ppm  i 

A.\ 

-41 1 5ST3.'  4 

-s.2  -'<iu 

AY 

-"1 

1 1.  s  -2s  .5 

AZ 

-91  r'T  '  YS 

-'3  -1" 

Slope  Distance 

i  M4  'i  i 

-<i  <\  ~  ~  2 

Table  21.  CASE  I:  DATA  SEGMENTS  BEFORE  AND  AFTER  REPROCESS¬ 
ING 


Dav  Segment 


Before 

After 

Rejected 

Obs 

Slope 
Distance 
( m) 

Rejected 

Obs 

Slope 
Distance 
t  m) 

Case  1 

14 

1 230050. 192 

0 

1230045.846 

17 

12302GS.9SS 

0 

1230045.294 

5 

1230047.276 

0 

1230047.020 

4i 

1230227.492 

0 

1230045.717 

44 

1230093.548 

0 

1230045.915 

6 

1230048.636 

0 

1230045.326 

23 

1230044.230 

12 

1 230046. 146 

* 

123O042.920 

0 

1230045.972 

Case  2 


m 

“S 

56 

1230074.209 

0 

1 230045. 2S2 

3.5  1 

Case  3 

_ 1 _ 

37, 

1230074.209 

. - 

_ a _ 

1230045.373 

3.5 

Case  4 

81 

1 

7S 

22 

1230034.198 

0 

1230045.082 

3.5 

SI 

2 

169 

47 

1230041.504 

11 

1230045.031 

3.5 

S2 

1 

62 

20 

1230039.213 

0 

1230045.024 

3.5 

Case  5 

S7 

2 

159 

43 

1 230043.41  S 

3 

1230045.177 

3.5 

OS 

l 

S6 

24 

1230043.047 

0 

1230045.472 

3.5 

14 

I 

52 

*  1 

1230051.620 

0 

123O04S.<m  (0 

99.0 

*"IT< 


currier  phase  time-tags.  All  the  Case  1  reprocessed  segments  were  improved  by  forcing 
the  processing  program  to  accept  more  carrier  phase  observations  by  increasing  the  edit 
multiplier. 

Because  of  the  success  of  improving  the  solutions  by  the  pseudorange  pre¬ 
processing.  all  segments  for  three  days  (  276,  28~.  and  294)  were  reprocessed  to  see  if  their 
already  acceptable  solutions  would  improve.  The  solutions  either  did  not  change  or 
became  slightly  less  accurate. 

F.  COMPARISON  WITH  PREVIOUS  STUDIES 

Remondi  [19S4)  concluded  that  100  ppm  accuracy  was  achievable  in  30  minutes  and 
10  ppm.  in  one  hour  -  regardless  of  baseline  length  -  based  upon  single-frequency,  single 
difference  solutions  using  precise  ephemerides  over  baselines  less  than  100  km.  To  study 
the  required  observation  times  for  various  accuracies,  he  partitioned  his  data  into 
15-minute  observation  spans.  While  decorrelated  triple  differences  would  be  capable  of 
performing  relative  geodesy  at  1  ppm  level,  correlated  triple  differences  may  provide  5-10 
ppm  or  better  [Remondi,  1985]. 

Some  long  baseline  surveys  include: 

Bock  et  al.  [19S4|  achieved  2  ppm  accuracy  in  closure  of  a  transcontinental  net  aided 
by  external  atomic  clocks  using  10  hours  of  single  frequency,  single  difference  obser¬ 
vations. 

Cannon  et  al.  [1985]  achieved  0.1  to  0.7  ppm  repeatability  using  single  frequency 
non-differenced  network  solutions  without  ionospheric  correction  for  l?00-km 
baselines  between  two  California  stations  and  Calgary.  Alberta.  Canada  over  a  four- 
day  period.  The  range  in  repeatability  was  attributed  to  satellite  clock  errors.  They 
also  found  that  the  largest  errors  occurred  in  the  AX  component. 

Mader  and  Abell  [1985]  found  that  single  frequency  long  baseline  GPS  results  using 
single  differenced  carrier  phase  observables  agreed  to  one  ppm  with  VLB1  measured 
baseline  lengths.  Their  2-day  repeatability  was  0.24  ppm. 

Goad  et  al  [1985]  measured  a  1302-km  baseline  between  California  and  Texas  to  better 
than  1.0  ppm  compared  to  the  VLB  I  measured  distance,  and  their  two-day  repeat¬ 
ability  was  better  than  0.5  ppm. 

Bertiger  and  Lichten  [1987]  achieved  6  parts  per  billion  repeatability  over  a  l?14-km 
east-west  baseline  over  separate  four  and  seven-day  periods  using  dual-frequency  re¬ 
ceivers  with  a  fiducial  network  for  orbit  determination,  and  some  water  vapor 
radiometers. 

Except  for  the  Bertiger  and  Lichten  [1987]  results,  the  Case  1  results  are  competitive 
with  the  above  studies  despite  the  inherent  low  resolution  and  the  added  burden  of  the 
correlations  of  the  correlated  triple  differences.  A  possibility  for  the  optimum  results  of 
Case  1  is  that  a  good  solution  is  insured  by  setting  both  receivers  to  their  known  WGS84 
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coordinates.  For  real-time  survevs,  one  end  of  the  baseline  will  usualiv  not 
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to  one  meter  accuracy. 


During  the  eight-week  observation  period,  the  third  4')0u$.\  was  delivered  to  NFS. 
The  antenna  was  installed  on  the  same  rooftop,  but  approximately  35  m  to  the  north 
of  the  \lonterev  antenna.  Two  days  of  data  usine  Monterey  as  the  reference  station 


were  used  to  fix  the  location  of  the  new  antenna.  On  Day  294.  the  new  receiver  s  coor¬ 


dinates  were  offset  3”  m.  equivalent  to  3d  ppm  baseline  error.  Die  difference  between 
the  olfset  antenna  and  Monterey  solutions  for  the  baseline  to  Sand  Point  is  o.“  ppm  for 
the  first  segment  <  F inure  1 1 1. 


Figure  1 1.  Results  of  station  2  offset 


For  the  remaining  segments  the  difference  is  about  0.05  ppm.  The  first  segment  will 
be  more  affected  by  receiver  offset  because  the  small  number  of  carrier  phase  observa¬ 
tions  will  be  sensitive  to  the  initial  position  estimate  provided  by  the  C  A  code.  The  C  A 
code,  in  turn,  will  be  sensitive  to  the  receiver  coordinates  because  the  receiver  will  com¬ 
pute  the  slant  range  to  the  satellites  using  those  coordinates. 


G.  MISSING  EPOCHS 


Missing  epochs  prevented  a  sixth  case  study  in  which  the  number  of  observations 
are  reduced  by  flagging  the  data.  The  sixth  case  study  was  tc  have  studied  the  effects 
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Figure  12.  Fraction  of  available  triple  difference  observations 

of  reducing  the  number  of  observations  without  changing  the  geometry.  Several  days 
of  data  were  processed  in  which  only  even.'  other  epoch  was  to  be  used  in  the  processing. 
An  hour  of  data  was  needed  to  obtain  enough  triple  difference  observations  to  consist¬ 
ently  have  less  than  10%  of  the  triple  difference  observations  rejected  and  slope  distance 
<7  s  less  than  10  m.  This  inability  to  form  a  reliable  solution  indicated  a  lack  of  carrier 
phase  observations,  so  the  number  of  observations  were  counted. 

The  mean  number  of  observations  (accumulative)  available  for  each  segment,  as  a 
fraction  of  the  maximum  number  of  triple  difference  observations  remained  steady  until 
120  minutes  when  it  increased  (Figure  12).  The  maximum  number  of  triple  difference 
observations  is:  \trD  =  {n,  —  !)(«,  —  1)  where: 

Mro  maximum  number  of  triple  difference  observations 
n,  number  of  observations  epochs 

n,  number  of  observed  satellites 

Solutions  using  only  10  minutes  of  data  at  a  time  were  examined  to  investigate  the 
rise  in  available  observations  after  120  minutes.  A  curious  phenomenon  was  discovered 
that  while  SV  6  was  being  tracked  at  Sand  Point  between  live  and  ten  of  even'  40  epochs 
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were  lost  (Figure  13).  As  soon  us  SV  6  was  no  longer  being  tracked  at  Sand  Point,  no 
complete  loss  of  epochs  was  observed  <  figure  14). 


START  OATE/TIME:  1987/10/21  16:43:07.  DAY  OF  YEAR  294  TOW  319360. 

STOP  DATE/TIME:  1987/10/21  16:52:60.  DAY  OF  YEAR  294  TCW  319980. 

DATA  AVAILABLE 
STATION:  SAND  POINT 


SV  6  I  . . 

SV  9  | . 

SV11  |  . 

SV12  I  . 

SV13  I  . 

STATION:  MONTEREY 

SV  6  I . 

SV  9  | . 

SV11  | . 

SV12  | . 

SV13  I . 


Figure  13.  Missing  epochs  while  tracking  SV  6:  Each  dot  represents  one  obser¬ 
vation:  each  column,  one  epoch.  Sand  Point  is  missing  epochs  4,  8,  14, 
.  .  20,  30.  34,  and  40.  From  Irim640  output. 


START  DATE/TIME:  1987/10/21  17:03:04.  DAY  OF  YEAR  294  TOW  320580. 

STOP  OATE/TIME:  1987/10/21  17:12:60.  DAY  OF  YEAR  294  TOW  321180. 

DATA  AVAILABLE 
STATION:  SAND  POINT 

SV  9  | . | 

SV11  | . | 

SV12  | . | 

SV13  I . | 

SV  6  I  | 

STATION:  MONTEREY 

SV  9  | . 1 

SV11  1 . | 

SV12  1 . | 

SV13  I . I 

SV  6  I . | 


Figure  14.  Continuous  tracking  without  SV  6.:  No  missing  epochs  at  Sand  Point. 


The  phenomenon  is  not  multipath  interference  because  it  affects  all  satellite  at  each 
epoch.  A  simultaneous  outage  would  suggest  interference  from  some  other  emitter  at 
or  near  1540  MHz.  if  the  outages  continued  throughout  the  observation  session  rather 
than  being  limited  to  when  SV  6  was  being  tracked.  No  physical  cause  for  the  outages 


is  readily  available,  and  the  data  have  been  forwarded  to  Trimble  Navigation  Ltd  for 
further  evaluation. 


H.  DILUTION  OF  PRECISION  AND  RANGE  ERRORS 

The  case  studies  showed  that  the  accuracy  and  precision  of  the  slope  distance  were 
improved  by  using  observations  associated  with  the  infinite  PDOP  peak  rather  than  ei¬ 
ther  the  PDOP  minimum  or  the  secondary  PDOP  rise  near  the  end  of  the  observing 
sessions,  and  that  improvements  in  the  slope  distance  accuracy  could  be  achieved  using 
fewer  observations  when  observations  are  made  through  the  infinite  PDOP  peak.  Pre¬ 
cision  and  accuracy  of  the  solutions  are  dependent  on  geometry  (fixed  with  respect  to 
Case  1  time)  and  error  sources  rather  than  the  number  of  observations  (fixed  with  re¬ 
spect  to  length  of  observing  period). 

The  mean  slope  distance  errors  were  governed  by  the  ionospheric  and  tropospheric 
errors  opposing  each  other.  Early  and  late  observations  had  positive  mean  errors  while 
the  errors  about  the  infinite  PDOP  peak  were  negative.  Tropospheric  errors  are  partially- 
corrected  by  the  modified  Hopfield  model  while  the  ionospheric  errors  remain  uncor¬ 
rected,  and  as  the  satellites  reach  the  infinite  PDOP  peak  the  full  impact  of  the 
ionospheric  errors  project  upon  the  baseline. 

The  role  reversal  of  the  PDOP  peak  from  point  positioning  to  relative  positioning 
requires  some  investigation.  Because  PDOP  reflects  the  strength  of  the  satellite  geom¬ 
etry  for  a  single  station,  the  investigation  should  begin  with  the  satellite-baseline  (two 
station)  geometry  by  examining  a  simple  projection  of  errors  upon  a  baseline  using  a 
single  satellite  whose  track  is  directly  over  the  baseline  (Figure  15). 

From  the  law  of  cosines: 

b\  =  p]  +  p\  ~  2p\Pi  cos  y  (6.1) 

where: 


bT  true  baseline 

p,  range  from  station  1  to  the  satellite 
p2  range  from  station  2  to  the  satellite 


Figure  15.  Relative  positioning  geometry:  Not  drawn  to  scale. 


Similarly  we  define  a  measured  baseline,  bM: 

bit  =  (Pi  +  £pi)2  +  (Pi  +  Epi )2  “  2(p,  +  ^,)(p2  +  Ef>2)  cos  y  (0.2) 


where: 

£,,  error  in  p, 
£,2  error  in  p, 


Subtracting  the  square  root  of  (6.1)  from  the  square  root  of  (6.2),  and  letting 
cosy  =  1,  gives: 


br  V  b\  b\ 


(6.3) 


where  AE,  =  Efi  -  Ef[  and  Ap  =  p}  -  p,  .  From  (6.3),  Ap  acts  as  an  amplification  factor 
for  the  difference  in  the  range  errors  (  A£,  ).  As  a  satellite  approaches  the  mid-point 
of  a  baseline,  Ap  -*  0  and  br>>  A E,. 


\  "w  -s.  s  *.  % 


Similar  uiguments  can  be  extended  to  point  positioning  living  two  vjtellnes  at  one 
epoch  or  two  epochs  of  one  satellite  oy  inverting  the  gcometrx  of  Figure  15.  but  the 
reduction  in  range  errors  would  apply  to  point  poMt.cnmg  as  well. 

1  lie  etiect'  ot  .  have  been  ignored  because  ,  is  small  1  rota  the  law  of  sines  : 

"  sin  H 


where  0  =  e 


e  is  the  elevation  unde  to  the  satellite  and  i  is  the  vertical  ande  be¬ 


tween  the  slope  d. stance  and  the  local  horizon.  I  xpanding  sin  t)  as 
sin  »’  cos  v  cos  *•  sin  v  .  and  :  sin  e  cos  •.)>>(  cos  e  sm  .  i  because  cos  .  ~  1  i  r  is  .s’  for 
the  lds'»  km  baseline  i.  then  sin  h  can  be  approximated  by  sm  e  Approximating  h  by 
1230  km  and  >>  In  2immn>  km.  then,  by  the  la  w  of  sines.  y  is  approximately  1°  when  »’  is 
15°;  by  the  law  of  cosines,  is  approximately  2.5°  when  the  satellite  is  at  the  baseline 
midpoint. 

With  cos  and  cos  2.5’iO  hhhl.  a  very  insignificant  change,  most  of  the 

deviation  of  the  measured  baseline  from  the  true  baseline  must  come  from  the  error  dif¬ 
ference  factor.  Most  importantly  having  satellites  at  high  elevation  angles  with  re¬ 
spect  to  the  termini  of  the  baseline  is  not  detrimental  to  relative  positioning  as  high 
elevation  angles  are  for  point  positioning. 

A  second  area  of  investigation  is  the  path  dependent  errors.  The  ionospheric  and 
tropospheric  errors,  approximated  by  cosecant  functions  of  the  elevation  angle,  are 
minimised  at  high  elevation  angles.  Simultaneous  high  elevation  angles  lor  both  stations 
is  another  by-product  of  high  PDOP. 

Finally,  the  effects  of  differencing  must  be  explored.  Qualitatively,  when  the  satel¬ 
lites  are  bunched  dose  together  at  high  PDOP  (Figure  3),  some  of  the  correlation  in 
errors  that  was  lost  in  lengthening  the  baseline  is  restored.  The  change  in  tropospheric 
and  ionospheric  errors  with  time  will  also  be  minimized  because: 


d  esc  e 
dt 


—  cos  e  d  e 
dt 


and  high  elevation  angles  minimize  Equation  (6.5).  Such  time  derivative  minimization 
of  errors,  will  greatly  affect  the  triple  difference  because  the  triple  difference  is  itself  a 
time  difference. 

The  sum  of  the  elevation  angles  for  ail  satellites  at  both  stations  (Table  22)  reaches 
a  maximum  after  70  minutes  of  observations,  which  is  ten  minutes  later  than  the  PDOP 


'  1  “Vv> 


j~L«  >J 


peak.  The  minimum  variance  is  achieved  in  the  SO  to  90  minute  Case  l  tune  The  run¬ 
ning  mean  of  the  sum  of  the  elevation  angle  peaks  at  the  inn  nunute  observing  tune 
while  the  running  mean  elevation  angle  readies  a  local  maximum  in  the  8<>  to  Inn  minute 
;ange.  I  he  running  means  provide  a  better  estimator  of  the  behavior  of  the  a^uracv 
and  preci'ion  because  of  the  accumulation  of  observations. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Following  a  recommendation  by  Remondi  [U)S4|  that  more  testing  ot  the  accuracies 
of  triple  dillerence  carrier  phase  measurements  was  needed.  1  studied  the  optimized  (/PS 
observation  times  required  to  achieve  geodetic  accuracies  by  partitioning  observation 
periods  for  a  1230-km  baseline  into  ten-nunute  segments  and  changing  the  length  and 
starting  time  of  observations.  Accuracy  was  determined  by  comparing  measured 
baseline  values  with  reference  v  alues  obtained  by  locating  the  ends  of  the  baselines  from 
very  precise  VLBI  horizontal  control  points  using  double  dillerence  GPS  carrier  phase 
measurements  over  short  baselines. 

1  he  1230-km  baseline  was  directly  measured  for  2S  days  over  a  period  of  eight  weeks 
using  Trimble  4<)tx.iSX  receivers.  The  Trimble  supplied  Tnnnec  software  was  used  to 
process  the  carrier  phase  measurements.  Using  broadcast  ephemcrides  in  a  lived  orbit, 
triple  dillerence  carrier  phase  solution  with  no  ionospheric  corrections  and  a 
tropospheric  delay  model  that  used  only  surface  meteorological  data,  I  At  part  per  million 
ippm)  accuracy  m  the  slope  distance  was  achieved  for  any  observing  period  with  a  day- 
to-day  repeatability  better  than  1.0  ppm  (2o).  The  AX.  AY,  and  AZ  (the  components 
ot  the  baseline  parallel  to  the  WGSS4  axes)  achieved  accuracies  better  than  10.0  ppm  lor 
any  observing  period.  AY  and  AZ  repeatabilities  weie  better  than  10.0  ppm  for  any 
observing  period,  while  AX  required  30  minutes  of  observations  to  reach  10.0  ppm. 

I  livid  not  expected  to  achieve  1.0  ppm  accuracy  m  the  slope  distance  measurements 
because  Remondi  [  1 9S4]  had  suggested  that  1.0  ppm  required  dual  frequency  measure¬ 
ments  using  highly  accurate  orbit  information  and  water  vapor  measuren  ents,  and  my 
solutions  are  further  burdened  by  the  correlations  of  the  triple  differences.  The  1.0  ppm 
accuracies  were  on  the  order  of  the  uncorrelated  single  frequency  results  of  previous  long 
baseline  GPS  survey  s:  Cannon  et  al.  [1985],  Bock  et  al.  [1984],  Goad  et  al.  [19S5],  and 
Mader  and  Abell  [1985]. 

The  unexpected  1.0  ppm  accuracy  can  be  attributed  to  low  ionospheric  activity  be¬ 
cause  of  the  orientation  of  the  baseline,  the  time  of  year,  and  the  minimum  in  the 
sunspot  cycle  and  solar  activity  during  the  observing  period.  Ephemerides  errors  are 
considered  low  because  of  the  low  age  of  the  ephemerides.  Those  ameliorating  factors 
must  be  considered  before  apply  ing  the  results  of  these  case  studies  to  planning  surveys 


on  different  baselines  at  different  times  ol  Jay  and  vear.  Ionospheric  errors  were  on  the 
order  ol  the  tropospheric  errors.  The  superposition  of  the  opposing  ionosphere  (short¬ 
ening!  and  tropospheric  (lengthening!  errors  reduced  mean  slope  distance  eriors. 

Slope  distance  errors  and  observation  periods  were  reduced  a  hen  Cl  PS  observations 
started  near  the  infinite,  symmetric  PDOP  peak.  I  he  reduction  of  etrors  about  the 
POOP  peak  can  be  attributed  to  the  smniltunecns  mmiira/ation  of  i.i.iee  ei:>'rs,  the 
projection  of  range  errors  onto  a  baseline,  satellite  and  epoch  separation.  A  high  sum 
of  the  satellite  elevation  angles  lor  a  single  station  can  aho  be  used  to  dcteimme  lav  ar¬ 
able  observing  times,  as  simultaneous  high  elevation  angles  correspond  to  hail  PDOP 
v  allies. 

My  results  reconfirm  Trimble  Navigation  s  | i^S7b)  recommendation  to  include  the 
infinite  PDOP  peak  in  the  observing  session;  however.  I  could  not  confirm  the  Federal 
Geodetic  Control  Committee  >  |  lVSoJ  proposal  to  stop  at  a  GDOP  of  5.0.  I  could  not 
identify  any  consistent  observing  stop  point  that  improved  accuracy  other  than  to  stop 
when  four  satellites  were  no  longer  available. 

When  less  titan  ?•>  minutes  of  data  are  used,  the  goodness  of  the  C  A  code  has  a 
great  died  upon  the  carrier  phase  solution  because  C  A  code  positions  are  used  to  es¬ 
timate  the  carrier  phase  solution  and  to  compute  the  tune  tags  of  the  carrier  phase 
measurements.  The  number  of  rejected  triple  dilference  measurements  provided  the  best 
indicator  of  the  quality  of  the  C  A  code  estimates  and  the  carrier  phase  solutions.  In¬ 
cidences  of  poor  accuracy  of  the  carrier  phase  solution  were  found  to  be  caused  by  poor 
C  A  code  estimates.  The  carrier  phase  solutions  could  then  be  improved  by  correcting 
the  C  A  code  solutions,  and  in  turn  the  carrier  phase  time  tags,  by  pre-processing  the 
C  A  code  measurements  or  accepting  more  carrier  phase  measurements.  Pre-processing 
the  C  A  code  measurements  did  not  improve  the  accuracy  of  the  carrier  phase  solutions 
when  the  number  of  rejected  triple  differences  was  less  than  10°  o  or  when  the  solution 
standard  deviation  of  the  slope  distance  was  less  than  10  m. 

The  long-term  average  of  the  best  pseudorange  solutions  produced  results  compa¬ 
rable  to  the  mean  carrier  phase  solution,  except  that  the  pseudorange  solutions  had 
v  ariances  several  times  the  magnitude  of  the  carrier  phase  solutions. 

B.  RECOMMENDATIONS 

The  results  of  these  case  studies  should  be  tested  on  surveys  conducted  with 
baselines  of  vary  ing  orientation  and  length  and  with  different  satellite  configurations. 
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Currently,  TRIMl'EC  dees  not  support  the  use  of  precise  ephemcrides  or  the  com¬ 
putation  ot  uncorrelated  triple  dill'erences.  Should  those  capabilities  be  implemented  m 
the  future,  or  provided  by  third-party  software,  the  data  should  be  reprocessed  and  an¬ 
alyzed  to  isolate  the  tropospheric  and  lonosphenc  error  contributions  and  to  improve 
the  baseline  component  results.  Both  of  those  capabilities  would  allow  better  insight 
into  the  olivets  of  geometry  and  the  tropospheric  and  ionospheric  errors. 

The  usual  role  of  the  correlated  triple  dillerence  is  to  aid  in  fixing  cycle  slips  when 
the  receiver  loses  lock  and  to  provide  the  initial  estimates  for  the  double  dillerence 
least-squares  processing.  The  current  data  should  be  reprocessed  in  the  double  differ¬ 
ence  mode  to  determine  whether  the  correlated  triple  difference  solution  provides  sulli- 
ciont  accuracy  to  fix  the  cycle  slips  and  estimate  the  initial  integer  ambiguities  in  long 
baseline  surveys. 

The  standard  values  used  to  estimate  accuracv  could  be  improved  by  a  more  rigor¬ 
ous  fixing  of  the  antenna  locations.  Both  sites  should  be  subjected  to  a  network  ad¬ 
justment.  either  locally  or  simultaneously  (possibly  in  conjunction  with  the  NGS  VLBI 
crustal  motion  studies).  The  extra  baselines  required  for  a  network  solution  could  be 
measured  using  the  third  Naval  Postgraduate  School  receiver. 

Studies  should  be  conducted  to  determine  the  effects  of  using  meteorological  obser¬ 
vations  far  removed  from  the  Sand  Point  antenna  site.  A  temporary  meteorological 
station  could  be  set  up  near  the  antenna  site  (possibly  in  cooperation  with  the  Weather 
Service  Forecast  Office  (WSFO)  located  on  the  grounds  of  the  Western  Regional  Center 
or  with  the  nearby  U  niversity  of  Washington).  Future  Naval  Postgraduate  School  GPS 
surveys  would  be  aided  by  portable  meteorological  instruments,  preferably  that  made 
digital  records  of  the  temperature,  humidity,  and  pressure. 

Because  of  the  importance  of  elevation  angles  to  accurate  results,  studies  should  be 
conducted  to  determine  the  performance  of  the  TRIMVEC  supplied  tropospheric  mod- 
els(  modified  Hopfield  and  Marini  )  at  various  elevation  angles. 

While  this  study  has  shown  the  utility  of  the  single-siaiion  PDOP.  a  more  complete 
DOP  may  provide  a  better  satellite  selection  aid.  Such  a  complete  DOP  (  or  Dilution 
Of  Relative  Position  (DORP))  should  incorporate  covariances  for  baseline  components, 
satellite  orbital  errors,  receiver  timing  errors,  ionospheric  and  tropospheric  delays,  un¬ 
certainties  in  the  reference  station  coordinates,  and  cross-covariances.  The  DORP 
would  also  be  formulated  with  respect  to  the  type  of  differencing  to  be  used. 
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Lastly,  the  amount  of  data  processed  for  these  case  studies  was  limited  by  the  speed 
of  the  microcomputer  and  the  necessity  of  transferring  the  results  to  the  Vn.d  If  st- 
graduate  School  mainframe  for  analysis.  As  more  OPS  surveys  will  he  conducted  in  the 
future,  the  demands  lor  processing  power  will  increase.  It  will  become  imperative  that 
the  processing  software  be  ported  to  the  mainframe  or  to  a  netwoiked  num-compu'.er. 
I  hose  computers  w.ll  allow  multiple  access  to  tiie  software  as  well  as  speeding  the 
processing.  As  an  interim  measure,  the  RAM  of  the  current  ensemble  of  microcomput¬ 
ers  should  be  increased  to  several  megabytes.  1  his  increase  in  RAM  will  allow  the  use 
of  RAM-Jisks  which  will  speed  up  the  data  loading  -  which  is  the  most  time  consuming 
of  the  processing  procedures. 


APPENDIX.  BATBLD.BAS  LISTING 

10  REM  PROGRAM  BATBLD. BAS  BY  R.  BOUCHARD  NOV87 

20  REM  PROGRAM  BATBLD:  BUILDS  THE  BATCH  FILE  FOR  TRIMVEC  PROCESSING. 

30  REM  READS  WX  OBS  FROM  SEA.  MET  FILE 

40  REM  DAY  283  REQUIRES  ITS  OWN  BATBLD  PROGRAM. 

50  INPUT  "batch  file  name";OFL$ 

60  SP$="  " 

70  S  l$="comniand  /c  tbf  nodd.  tem 
80  OPEN  OFL$  FOR  APPEND  AS  02 

90  DIM  A$(  18),JD$(29)  ,MM$(  29)  ,DD$(  29)  ,PI!(  29)  ,PM(  29) 

95  REM  initialize  Julian  Day  array 

100  JD$( n="272":  JD$( 2)="273":  JD$( 3)="274":  JD$(4)="275" 

:  JD$( 5 )="276" 


110  JDS(6)="279":  JD$( 7 )="280":  JD$( 8)="28l":  JD$(9)="282" 

:  JD$( 10)="287" 

120  JD$( 1 1 )="290":  JD$( 12)="293":  JD§( 13)="294":  JD$( 14)="295" 


JD$(  15)="297" 


130 

JD$( 16)=" 

300": 

JD$( 17 )="30 l": 

JD$( 18)="302":  JD$( 19)="307" 

140 

JD$( 20 )=" 

308": 

JD$(21)="309": 

JD$( 22 )=" 3 14" 

145 

JD$(23)=" 

315": 

JD$(24)="317": 

JD$( 25 )=" 318" 

147 

J0$( 26 )=" 

325": 

JD$( 27 )="329" 

:  JD$(28)="321" 

150  MM$(1)="09":  MM$(2)="09" 

160  FOR  1=3  TO  18:  MM$(I)="l0":  NEXT  I 
170  FOR  1=19  TO  28:  MM$(I)=”ll":  NEXT  I 

180  DD$( n="29":  DD$(2)="30":  DD$(3)*"0l":  DD$(4)="02":  DD$(5)="03" 
:  DD$(6)='06":  DD$( 

7)="07":  DD$(8)="08":  DD$(9)="09":  DD$( 10)=" 14":  DD$(ll)="l7" 

:  DD$(J2)="20" 


A$(  2) 
A$(  8) 


»"a08": 


A$(9)="a09":  A$( 10)="al0":  A$(ll)="all" 
:  A$( 15)="al5":  A$( 16)="al6" 


190  DD$(13)="2l":  DD$( 14)="22":  DD$( 15)="24":  DD$( 16)="27":  DD$(17)="28" 
200  DD$( 18)="29":  DD$( 19)="03":  DD$( 20)="04":  DD$( 21)="05":  DD$(22)="l0" 
205  DD$( 23)=" 1 l":  DD$( 24)=" 13":  DD$( 25 )=" 14" 

207  DD$(26)="21":  DD$( 27 )="25":  DD5(28)="l7" 

210  A$( 1 )="a0 l":  A$(2)="a02":  A$(3)="a03":  A$(4)="a04":  A$(5)="a05" 

:  A$( 6)="a06" 

220  A$( 7 )="a07":  A$(8)="a08":  A$(9)="a09":  A$( 10)="al0":  A$(ll)="all" 

:  A$(12)="al2" 

230  A$(l3)="al3":  A$( 14)="al4":  A$( 15)="al5":  A$(16)="al6" 

:  A$( 17)="al7":  A$(18)="al8" 

240  FOR  1=1  TO  7:  PH(I)=19:  NEXT  I 
250  FOR  1=8  TO  13:  PH(I)=18:  NEXT  I 
260  FOR  1=14  TO  21:  PH(I)=17:  NEXT  I 

270  FOR  1=22  TO  25:  PH(I)=16:  NEXT  I 

275  FOR  1=26  TO  28:  PH(I)=15:  NEXT  I 

280  PM(1)=34:  PM(2)=31:  PM(3)=26:  PM(4)=23:  PM(5)=18:  PM(6)=5:  PM(7)=1 
290  PM( 8)=58:  PM(9)=54:  PM(10)=33:  PM(11)=22:  PM( 12)=8:  PM(13)=3 
:  PM( 14)=59 

300  PM( 15)=51:  PM( 16)=39:  PM(17)=35:  PM( 18)=31:  PM( 19)=14 
310  PM( 20)=6:  PM(21)=2:  PM(22)=42 

315  PM( 2 3 ) =3 7 :  PM(24)=29:  PM(25)=24:  PM(26)=56:  PM(27)=39:  PM(28)=8 
317  REM  begin  building  the  batch  file  for  each  day. 

320  FOR  IK=1  TO  28 
330  OPEN  "i", 1, "sea. met" 


260  FOR  1=14  TO  21 
270  FOR  1=22  TO  25 
275  FOR  1=26  TO  28 


54 


)  IN'F$="sa"+JD$f IK)+"  ma"+JD$(IK) 

)  S2$=S 1S+INFS 
)  PRINT  "FCR  JD:  JDS( IK) 

)  INPUT  ''start  hour";  SH 
)  INPUT  "start  minute"; SM 
)  PRINT  JDS ,PH( IK) ,PM( IK) 

)  SilS=STRS(. SH) 

)  SM$=STR$(SM) 

)  FC=SM/bO 
)  1  =  1 

.  REM  Find  the  WX  OB 

i  INPUT,?  1,  IDS.IHS  ,P1(  I)  ,T1(  I)  ,R1(  I)  ,P2(.  I)  ,T2(  I)  ,R2(  I) 
i  IF( f IDS=JD$( IK) )AND( SH=VAL( IHS j ) )  THEN  470 
i  GOTO  440 

'  Pl( I )=P1( I ) -2  :  T2(I)=(T2(I)-32)*5/9 

i  1=1+1 

INPUTS  1 , IDS , IHS , Pl( I ) ,T1( I) ,R1(I),P2(I),T2(I),R2(I) 
I F C IDS <> JD$ ( IK))  THEN  540 
REM  Subtract  2  mb  from  SEA  ob 
P I ( I )=P1( I ) -2! 

REM  Convert  NPS  Temp  to  Celsius 
T2(I)»(T2(I)-32)*5/9 
GOTO  480 
IE=I - 1 

REM  interpolate  WX  OBS  TO  GPS  START  TIME 

PR1(1)=  Pl( 1)+(P1(2)-P1( 1))*PC 

PP-2(  1)*P2(  1)+(P2(2)-P2(  1))*PC 

TR 1  ( 1 )  =T  1  ( 1)+(T1(2)-T1(1))*PC 

TR2(  1 )=T2(  1)+(T1( 2) -Tl( 1) )*PC 

RP 1C 1)=R1(  1)+(R1(2)-R1( 1) )*PC 

RR2(  1)=R2( 1 )+( R2( 2) -R2( 1) )*PC 

REM  Compute  running  mean  of  first  hour 

PR1(2)=  C  PR  1 ( I )+Pl( 2) )/2 

PR2(2)=(PR2( l)+P2(2))/2 

TR 1 ( 2 ) =( TR 1 ( l)+Tl(2))/2 

TR2( 2)=(TR2( 1)+T2( 2) )/2 

RR1(2)=(RR1( l)+Rl(2))/2 

RR2( 2)=( RR2( 1)+R2( 2) )/2 

WT=1-PC 

REM  Compute  remaining  running  means 

FOR  1=  3  TO  IE 

IA-I-i 

IWT=1+WT 


720  PR1(I)=(PRI( IA)*WT  +  ( ( Pl( IA)+P1( I ) )/2) )/( IWT) 
730  PR2( I )=( PR2( IA)*WT  +  ( ( P2( IA)+P2C I ) )/2) )/( IWT) 
740  TR1( I )®(TR1( IA)*WT  +  ((Tl(IA)+Tl(I))/2))/( IWT) 
750  TR2(I)=(TR2(IA)*WT  +  ( (Tl( IA)+T2( I) )/2) )/( IWT) 
760  RR1(I)=(RR1(IA)*WT  +  ((Rl(IA)+Rl(I))/2))/( IWT) 
770  RR2( I )— ( RR2C IA)*WT  +  ( (R2( IA)+R2( I ) )/2) )/( IWT) 
730  WT=IWT 
790  NEXT  I 
800  M1=SM 
810  H=SH 
820  C=0 
830  1=2 

840  DHR=PH( IK) -SH:  DMIN=PM( IK) -SM 


55 


850  IF(DMIN<0)  THEN  DHR=DHR- 1:  DMIN=DMIN+60 
860  NP=DHR*6+(DMIN/  10) 

865  REM  Interpolate  hourly  running  means  to  10  minute  intervals 
870  FOR  I L=  1  TO  NP 
880  M1=M1+10 

850  IF(M1>60)  THEN  Ml=Ml-60:  H=H+1:  1=1+1 
500  WT=Ml/60 

910  IF(M1=60)  THEN  Ml=Ml-60:  H=H+1:  1=1+1 

920  NP=M 1 

930  HP=H 

990  IWT=1-WT 

950  IA= I -  1 

960  PV1=PR1( IA)*( I WT ) +PR 1(1) *WT 
970  Pl,T=INT( ( 10*PV1)+. 5)/10 
980  PW2=PR2( IA)*( IWT)+PR2( I)*WT 
990  PW2=INT( ( 10*PW2)+. 5)/ 10 
1000  TW1=TR1( IA)*( IWT)+TR1( I)*WT 
1010  TW1*INT((  10*TW1)+.  5)/10 
1020  TV2=TR2(  IA)'V(  IWT)+TR2(  I  )*WT 
1030  TW2=INT( ( 10*TW2)+. 5)/ 10 
1040  RW1=RR1( IA)*( IWT)+RR1( I)*WT 
1050  RW1=INT((10*RW1)+. 5 ) / 10 
1060  RW2=RR2( IA)*( IWT)+RR2( I )*WT 
1070  RW2=INT(( 10+RW2)+. 5 ) / 10 
1080  REM  convert  to  strings 
1090  PS1$=LEFT$(STR$(PW1) ,7) 

1100  PS2$=LEFT$(STR$(PW2) ,7) 

1110  TS1$=L£FT$(STR$(TW1) ,5) 

1120  TS2$=LEFT$( STR$(TW2) ,5) 

1130  RS1$=LEFT$(STR$(RW1) ,5) 

1140  RS2$=LEFT$( STR$( RW2) , 5 ) 

1150  HS$=STR$(HP) 

1160  S3$=S2$+SP$+A$( IL)+".  "+JD$(IK) 

1170  MS$=STR$( MP) 

1180  PRINT  S3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+SP$+MM$( IK)+SP$ 
+DD$( IK)+HS$+MS$+SH$+SM$ 

1190  PRINT^2,  S3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+SP$+MM$( IK) 
+SP$+DD$( IK)+HS$+MS$+SH$+SM$ 

1200  LPRINT  S3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+SP$+MM$( IK)+SP$ 
+DDS( IK)+HS$+MS$+SH$+SM§ 

1210  NEXT  IL 
1220  CLOSE  U 
1230  NEXT  IK 
1240  CLOSE  ‘41 
1250  END 
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